Jet quenching in hot strongly coupled gauge theories revisited: 
3-point correlators with gauge-gravity duality 



Peter Arnold and Diana Vaman 

Department of Physics, University of Virginia, 
Box 400714, Charlottesville, Virginia 22904, USA 
(Dated: October 26, 2010) 

Abstract 

Previous studies of high-energy jet stopping in strongly-coupled plasmas have lacked a clear 
gauge-theory specification of the initial state. We show how to set up a well-defined gauge theory 
problem to study jet stopping in pure AA=4 super Yang Mills theory (somewhat analogous to 
Hofman and Maldacena's studies at zero temperature) and solve it by using gauge-gravity duality 
for real-time, finite-temperature 3-point correlators. Previous studies have found that the stopping 
distance scales with energy as E 1 ^ (with disagreement on the gauge coupling dependence). We do 
find that none of the jet survives beyond this scale, but we find that almost all of our jet stops at a 
parametrically smaller scale proportional to (EL) 1 ^, where L is the size of the space-time region 
where the jet is initially created. 



I. INTRODUCTION AND RESULTS 



How far does a localized, high-energy excitation travel in a quark-gluon plasma before 
slowing down, stopping, and thermalizing? This question is of phenomenological interest 
for heavy ion collisions, and it has long been a problem of theoretical interest to calculate 
the result in various idealized situations. In a weakly-coupled gauge theory with massless 
partons, the stopping distance scales with energy as up to powers of logarithms, 

where E is the initial energy of a high-energy parton. 1 At the other extreme, investigations 
0-0] of strongly-coupled, large-A c , supersymmetric versions of QCD, using gauge-gravity 
duality 0413, have indicated that the maximum stopping distance scales like E 1 ^ 3 . In this 
paper, we revisit this gauge-gravity duality result. In particular, previous calculations have 
always specified the high-energy initial state using the gravity description (and in some 
cases relegated other matters of interpretation to the gravity description): there has not 
been a complete specification of the problem, from beginning to end, solely in terms of 
4-dimensional gauge theory. We will investigate what happens if a localized, high-energy 
excitation is created in the gauge theory, and the response later measured in the gauge theory. 
For our method of creating the initial excitation, we find that there is an additional scale 
characterizing the response: almost all of the excitation's conserved charge is deposited at a 
distance that scales with energy as E 1 ^ rather than E 1 / 3 and is sensitive to the initial spatial 
size of the excitation. Nonetheless, we will still see E 1 ^ 3 emerge as the furthest distance that 
any non-negligible fraction of the charge propagates before stopping and thermalizing. 

To be more concrete, we need to explain more precisely what we calculate. 



A. The problem 

The specific theory we study is pure A/"=4 super Yang Mills theory in the large iV c and 
large A = gy M N c limit. Readers needing a general introduction to the use of gauge-gravity 
duality to study finite-temperature physics in this and similar strongly-coupled theories 



should refer to Ref. 12 



We will follow the general philosophy of Chesler et al. 0, [ll| that the way to study stop- 
ping distances is to locally create a high-energy excitation and then measure the subsequent 
evolution of conserved charge densities such as energy or momentum density. A cartoon 
of this evolution is shown in fig. [TJ In contrast to previous studies at finite temperature 
T, we will give an explicit gauge theory prescription for creating the initial excitation. We 
excite the gauge theory plasma by turning on external sources, localized in space-time, that 
produce a high-energy state with nearly-definite energy and momentum. One could in prin- 
ciple use most any type of source that has a simple translation to the gravity dual theory 
under the AdS/CFT correspondence, but in this paper we will focus on an example where 
the sources couple to the global R-charge currents of the gauge theory. It will also simplify 
the analysis to use a source that is translation invariant in the two spatial dimensions trans- 
verse to the motion of the excitation, but localized in time and the third space direction. 



1 A specific calculation for QCD of the stopping distance at weak coupling in the high-energy limit may 
be found in Rcf. 1]. However, the scaling of this result was implicit in the early pioneering work of Refs. 
0, 0] on bremsstrahlung and energy loss rates in QCD plasmas. 
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Specifically, we modify the 4-dimensional field theory Lagrangian by 



C^C+f^, (1.1) 

where j" are the SU(4) R-charge currents of the theory and A c \ is a classical external source. 

We choose the external source to have the form of (i) a high-energy plane wave e lk ' x times 
(ii) a smooth, slowly varying, real- valued envelope function A L (x) localizing the source to a 
space-time region of size L. Specifically, 



where 



— e fk - x + h.c. 
L 2 



A L (x), (1.2) 



= (E, 0, 0, E) (1.3) 

is a very large light-like 4-momentum with frequency E ^> T; Ma is an arbitrarily small 
source amplitude; e is a transverse linear polarization, such as 

^ = (0,1,0,0); (1.4) 

and t 1 are Pauli matrices for any SU(2) subgroup of the SU(4) R-symmetry, with 

T ± = r 1 ±ir 2 . (1.5) 

(The motivation for the r + factor will be discussed below). A simple example of an appro- 
priate envelope function would be 

A L ( X ) = e ~U-o/Lf e -U^/L)\ ( L6 ) 

L should be chosen large compared to 1/E, so that the momentum components in the source 
are all close to (11. 3p . but small compared to the large stopping distance that we wish to 
study. To avoid the mental clutter of an over- abundance of scales, it is convenient (but not 
necessary) to consider L to be of order 1/T in what follows. 

The source ( II. 2p creates an excitation that carries energy, momentum, and R charge. 
We could subsequently track the densities of any of these conserved charges to study the 
evolution of the excitation. In this paper, we have chosen to study the evolution of the R 
charge density, specifically the large-time behavior (t ^> both T _1 and L) of 

</ 3) V)) Aci (1-7) 

if the system starts in thermal equilibrium at t = — oo. Here, the superscript "(3)" indicates 
the R charge current associated with r 3 /2 in the SU(2) subgroup referenced by (ll.5j) . and 
the subscript U A C \ indicates that the expectation is taken with the source term (II. ip present 
in the Lagrangian. Because of the r + in ( 11. 2p . the situation is analogous to the interaction 
of a quark-gluon plasma with an external, high-energy W + boson in a wavepacket of size 
L (and with a decay time of order L), as depicted in fig. [2j 2 The W + boson will leave 
behind an excitation that carries electric charge and the third component of isospin r 3 /2. 
Subsequently measuring the latter is analogous to (JTTTJ) . 



2 One could imagine strengthening this analogy by gauging an SU(2) subgroup of the 7V=4 super Yang 
Mills R symmetry with a very weak coupling constant j w < 1. The full SU(4) R symmetry is anomalous 
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FIG. 1: The space-time development of a conserved charge density carried by an initial high- 
energy excitation that is moving along the light cone and interacting with the thermal medium. 
The development transitions between a ballistic trajectory at early times to diffusion at late times, 
(a) shows sketches of density vs. x at a sequence of larger and larger times; (b) depicts the space- 
time region where the density is non-negligible. The red circle at the origin of space-time denotes 
the region of size L where the source (|1.2p that creates the initial excitation is non-negligible. 




FIG. 2: A very high energy W + boson decaying inside a standard-model quark-gluon plasma and 
producing high-energy partons moving to the right with net 3rd component of isospin, r 3 /2. In 
the context of AA=4 super Yang Mills, the u and d above represent adjoint-color fermions or scalars 
carrying R charge and, for strong coupling, should not be pictured perturbatively as in this picture. 



and so cannot be consistently gauged unless one adds yet other fields to the model to cancel the anomaly. 
But an SU(2) subgroup is not anomalous and could be gauged, provided one defines the currents of that 
subgroup appropriately. The currents that the gauge bosons would have to couple to would be slightly 
different than the usual currents defined in the AdS/CFT correspondence with holographic regularization 
because the latter treat all the R currents on an equal footing. This difference in currents reflects the 
difference between the covariant and consistent anomalies [3| ■ None of these distinctions actually matter 
in the current problem with source (jl.2l) . but we will simply avoid gauging any of the R currents so that 
we do not need to ponder these issues. 
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B. The result 



We will take the source amplitude Ma to be arbitrarily small so that we can treat the 
source term (11.11) as a small perturbation in our later analysis. A small-amplitude source 
will most of the time have no effect at all on the system, producing no excitation and no 
R charge. We can normalize away this case simply by dividing the average charge density 
distribution (j( 3 '°{x)) A by the average total charge produced by the source, 



= J d 3 x <j (3)0 (x)X 



2°>L 



In the case of a transverse-translational invariant source, such as will be studied in this paper, 
it is more appropriate to consider the charge per unit transverse area = /V±. 

The simplest way to express our final result is to give a charge-deposition function 0(x) 
which represents how much thermalized charge the high energy excitation leaves behind at 
each space-time point x. More specifically, Q(x) is the source term for the diffusion equation 
for the subsequent evolution of that charge, so that the late-time charge density is given by 3 



(d t -DV 2 ){j®°(x)) A ~Q™e(x), 



(3) 



;i.9) 



where the ~ here indicates that we are only resolving structure on distance and times scales 
large compared to the thermal wavelength ~ 1/T. For a strongly-coupled plasma, that is 
the hydrodynamic limit — the limit where the diffusion equation is applicable. The value of 
the R-charge diffusion constant D is 14 



D 



1 



2nT 



1.10) 



Our result is that, if one contents oneself with only resolving details on distance scales 
large compared to both the source size L and the thermal wavelength 1/T, then 



Q(x) ~2 5 L (x~ 




(ic 4 EL) 2 



c 4 EL 



(2ttT) s (x+)V * ^ (2ttTx+) 4 



E 1 / 3 



X 



X 



< J E 1 / 3 /(27rT) 4 / 3 ; 
>£ 1 / 3 /(2ttT) 4 / 3 . 

where 8(x + ) is the step function; x^ = x 3 ± x°; and ^(y) is a source-dependent function 
that suppresses \y\ 1, causing suppression of x + <C (EL) 1 / 4 / \2ixT) above. In the case of 
the Gaussian source f ll.6p . 

9{y) = e-^\ (1.12) 

The subscript L on 5l(x~) indicates that 5l(x~) is only an approximate delta function, with 
a width of order L. Approximating it as a true delta function is good enough if we are only 
interested in the hydrodynamic response on scales large compared to L. In (II. lip , the c's 
are numerical constants, given by 



r 2 (|) 

(2?r) 1 /2 



ci ~ 0.927, c 2 ~ 3.2. 



(1.13) 



3 The idea of defining and investigating 9 (x) has been taken from Chesler, Jensen, and Karch . However, 
in that paper, they applied it only to scales large compared to the stopping distance of the high-energy 
excitations. Here, we resolve all scales where hydrodynamics is applicable. 
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t 




(EL) V4 T 1 



(a) 



exponential 
fall-off 



E 1/3 T -4/3 



(EL) V4 T 1 

(b) 



£ 1/3 r -4/3 



FIG. 3: The deposition of charge, (a) shows the coefficient of the 5l(x~) in (II. lip for as 
a function of x + . (b) depicts the space-time points where O(x) is not exponentially suppressed, 
with lighter and lighter shading representing the algebraic (x + ) -9 fall-off of the strength of Q(x). 
In both figures, the axis are not represented linearly. 



We shall see later that C\ and C2 are determined by the first quasi-normal mode in the gravity 
description. 

A qualitative summary of the the x + dependence of Q(x) is shown in fig. |3] The deposition 
of charge furthest from the origin that is not exponentially suppressed is at 



Mi 



£1/3 
£4/3 ' 



;i.i4) 



which scales with energy as E 1 / 3 like the various results of Refs. |4j-l6|. However, at least on 
average, only a tiny 0([L/ (x 3 ) max ] 2 ) fraction of the total charge is deposited at this distance 
if we keep the source size L small compared to (x 3 ) max itself. Most of the charge is deposited 
at the much smaller distance scale 



(*s) 



dominant 



(EL) 



1/4 



T 



;i.i5) 



There is an important difference between (j!.14j) and the subset 0, 0] of earlier results 
which have studied jet stopping by studying the dynamics of classical strings on the gravity 
side. Ref. jsj (see also added massless J\f—2 fundamental- charge matter to supersym- 
metric Yang Mills and argued that one could study the stopping of excitations carrying 
the analog of baryon number by studying the stopping of moving classical strings in the 
gravity dual. Ref. 0] studied the pure jV=4 Yang-Mills theory and modeled gluon jets by 
the evolution of folded pieces of string in the gravity dual. In both references, the stopping 
distance was found to be of order 



(*s) 



E, 



string 



1/3 
string 



^1/6^4/3 



(1.16) 
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which is parametrically smaller than our corresponding scale (I1.14p by a factor of A~ 1//6 in the 
strong-coupling limit A — > oo. Formally, the origin of the factor of A _1//6 in the calculations 
based on classical strings is the fact that the string tension, and therefore the energy of 
the state represented by the string, is proportional to a/A, so that E^ mg oc A 1 / 6 . In our 
calculation, in contrast, we will not consider classical strings at all. We will just use the 
most basic, original elements of the AdS/CFT dictionary for relating gauge theory operators 
to classical boundary sources in the gravitational dual. 

Our maximal scale (11. 14j) does agree (including the absence of A) with the scale previously 
found by Hatta, Iancu, and Mueller jBJ, who specifically studied R charge excitations like 
we do. Gauge-gravity duality relates the SU(4) R currents to classical 5-dimensional SU(4) 
gauge fields on the gravity side. They studied how a wave solution of the gravity-theory fields 
would fall into the black brane horizon, and then they used rough, qualitative arguments 
to relate this behavior back to what happens in the 4-dimensional gauge theory. In this 
paper, we precisely relate the field theory problem we have outlined to a calculation in 
the gravity dual. Along the way of solving the gravity dual problem, we will eventually 
encounter the same sort of problem studied by Hatta, Iancu, and Mueller. However, we 
will find another scale to the problem, which was missed in their qualitative interpretation, 
and which corresponds to the scale (I1.15P at which almost all of the charge is deposited. In 
fact, we will see that the appearance of this scale from the gravity calculation is intimately 
related to the conservation of R charge in the 4-dimensional gauge theory problem after the 
source turns off (t 3> L). 

Before moving on, we record the generalization of our result (11. lip to the case of a generic 
source envelope, which is that 

Jdq_ |Aj 2) (9 +! g- 
4L*fdq + dq-0(-q + )\q + \\A™{q + ,q. 

where 9 is the step function and A^ 2 ^ is the two-dimensional Fourier transform 

Af{q + ,q_) = J ^^A L (x)e-^ ++ ^-\ (1.18) 

Provided the envelope function Al(x) is smooth on the scale L and falls sufficiently rapidly 
for ^> L, the qualitative conclusions are the same as for the Gaussian envelope. 



*(W0 = \.r L T,7" mi rrr. d-"> 



C. Jets in strong coupling 

We loosely use the term "jet" to refer to a spatially localized, high energy excitation 
that initially moves at nearly the speed of light through the plasma. There is a potential 
for confusion because sometimes people loosely summarize the (zero temperature) results of 
Hofman and Maldacena [THj] as indicating that there are no jets in strongly-coupled Af=4 
super Yang Mills. Hofman and Maldacena considered the case of an isotropic source localized 
in all four space-time dimensions, with 4-momentum narrowly peaked around k = (E, 0, 0, 0). 
In weak coupling, this source would predominantly produce two back-to-back partons, flying 
in opposite directions, and so produce a non-trivial angular distribution in the late-time 
energy-energy correlation function far away from the source. In the strong coupling limit, 
they found the opposite: there was no such angular correlation. Their source did not create 
a jet-like structure. 
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However, at zero temperature, one may change the appearance of something simply 
by boosting to a different reference frame. An expanding spherical shell of energy in the 
original frame looks like a single, slowly spreading, localized jet if one boosts by a very large 
amount in the x% direction. Effectively, this is what our source ( II. 2p does when we choose 
k = (E, 0, 0, E) in (II. 3p . The envelope function Al(x) causes a narrow spread in momentum 
k around (E, 0,0, E), and we will see later that it is the time-like subset {k^k^ < 0) of 
these momenta that produce our result (see also [5j). So the physics of the creation of 
our initial "jet" state is essentially an extremely boosted version of the zero-temperature 
problem studied by Hofman and Maldacena. 

Hofman and Maldacena related the measurement of one-point correlations (S(x)) of en- 
ergy density to the calculation of three-point correlations (0|O^£(a:)Ofc|0), where Oj, was 
the operator that created their initial state. In order to measure charge densities such as 
( II. 7p in this paper, we will similarly investigate three-point correlators (in our case at finite 
temperature) between (i) the measured charge density and (ii) operators associated with the 
creation of the source. However, for reasons of calculational simplicity that we will explain 
later, we have set up the problem in a slightly different way than Hofman and Maldacena 
did and so evaluate a slightly different type of 3-point correlator ordering. Specifically, we 
cast the problem in terms of retarded 3-point correlators. 



D. What follows 

In the next section, we set up the basic integrals that we will have to evaluate to obtain the 
charge density in terms of bulk-to-boundary propagators in AdSs-Schwarzschild 

space. Then we warm up to the task of evaluating these integrals in section UTTl by applying 
our method to the simpler case of zero temperature. In particular, we will make comparison 
with Hatta, Iancu, and Mueller [5( in section IIII El We move on to the finite-temperature 
case in section HVt where we derive our final result (II. lip . Finally, we conclude with some 
suggestions for future work in section [V] 



II. GENERAL SET-UP 



A. Notational Preliminaries 



We will use the form of the AdSs-Schwarzschild metric given by 



ds 2 = ^ 
4 

4 



-(-fdt 2 + dx 2 ) + ^- 
\2nTf 



u 



-f dt z + dx z ) + 



du 2 
1 



u 2 f 



du 



(2.1) 



where 



/ = 1 - (2ttT) 4 m 2 = 1 - u 2 , (2.2) 

R is the radius of the AdS space-time (which results will not depend upon), u=0 corresponds 
to the 4-dimensional boundary, and u—1 corresponds to the horizon. For space-time indices, 
we will use capital roman letters (/, J, ...) for indices in 5-dimensional space-time and Greek 
letters (fi, u, ...) for four-dimensional space-time. When we write a lower Greek index on a 



S 



4-momentum Q or polarization s, we will always mean that the index is lowered with the 
4-dimensional metric rj^ and not the 5-dimensional metric gu\ so 



Q/j, — i]/j.vQ i — f)ia>^ i (2-3) 

where r\ = diag(— 1, 1, 1, 1). 

Our conventions for light cone coordinates will be 

x — x —I — x . x _[_ — 2" x — ~~ . x — — 2" x — ~ - ^2.4^ 
and similarly 

? ± = 9 S ±9°, 9+ = !?- = *-=^, q- = lq + = q -^, (2.5) 

and so q^ = q + x + + q^x~ + q±- x±. When integrating over 4- momenta q, we will use the 
short-hand notation 

f f dAq f 2 dq+ dq ~ d2q± (0 a\ 

J q "'-JW~-J (2tt)« {2 - b} 

Where there is no opportunity for confusion, we will abbreviate (j^ 3 ^(x)) A i as 
and sometimes even as j^(x). 



B. Field Theory: 3-point functions 

To relate the response of 1-point functions such as energy density or R charge density 
^'( 3 )°(x)) to equilibrium n-point correlation functions, for small source amplitudes Ma-, one 
follows the same steps as in derivations of the fluctuation-dissipation theorem or Kubo 
formulas. Since it's relatively simple, we'll take a moment to review it here. Write the 
Hamiltonian as H(t) = H + 8H(t), where SH(t) are the small-amplitude source terms and 
H is everything else in the full Hamiltonian of the theory. If the system initially starts in 
equilibrium, before the sources turn on, then the later evolution of an observable O is given 
by 

(0(t)) H = Z, 1 tr (e- fiHo [U(t, -00)]* OU(t, -00)) , (2.7) 

where 

U(t, t ) = T exp ^-i jf dt' H(t'^J (2.8) 

is the evolution operator under H, and T is time ordering. Working in the interaction 
picture and expanding in powers of the small 5H, one finds 

(0{t)) H - (O) H0 = J dt, G R (t i; t) + ^J dh dt 2 G R (t u t 2 ; t) + ■ ■ ■ (2.9) 
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where the various Gr, are the equilibrium n-point retarded correlation functions, given in 
this case by 4 



iG R {h; t) = 6(t - h) ([0{t), 6H(ti)]) Ho , (2.10) 
i 2 G R (h,t 2 ; t) = 6(t - t 2 ) 6(t 2 - t x ) ([[0(t), 5H(t 2 )} 1 5H(t 1 )}) Hq 

+ 9(t - h) dih - t 2 ) ([[0(t), 6H(tx)], $H(h)]) Ho , (2.11) 

etc. In our problem, the first non-vanishing term is the one involving the 3-point correlator 
because the operator [j^°(x), j a ■ A^(xi)] has non-zero R-charge, causing the 2-point corre- 
lator Gr(ii; x) to vanish. However, the conclusion is more general than our specific example 
involving R charge currents. To create a high-energy excitation, the source 5H should have 
large momentum k. To measure later hydrodynamic behavior after the jet stops in the 
medium, one wants to examine relatively /oty-wavenumber components q of the late-time 
diffusing density 0(x). Because of this momentum mismatch between source and observ- 
able, the equilibrium two-point function Gn(k; q) will vanish by momentum conservation. It 
is only when we get to the three-point function that we first find a non- vanishing result. 5 
In our problem, (12. 9 j) manifests in detail as 

0' (3)M (*)>a c1 = I / *x x d A x 2 G^ b3)a ^(x 1 ,x 2 ;x)A a a Jx 1 )A% cl (x 2 ) 

= \ [ Gt m)a ^{Qx, Q 2 ; Q) A^iQx) A% cl {Q 2 ) e^{2n)H^{Q 1 + Q 2 + Q) 
JQ1Q2Q 

(2.12) 

in the limit of arbitrarily small source amplitude Ma- 

Now comes a crucial argument that we will use repeatedly. In studying the response 
(I2.12p to determine how far the jet travels, we will not care about the detailed structure on 
small distance and time scales such as 1/E or even 1/T. It would be perfectly adequate to 
look at a smeared response such as 

, C p~ (^ x °) 2 /^smcarp — I 1 2 /^smcar 

0' (3)At (-))_ d - / (jW»(x + Ax)) Aci , (2.13) 

j ^ smear 

where the smearing distance 4mear is chosen large compared to microscopic scales such as 1/E 
and 1/T but small compared to scales we're interested in resolving, such as the stopping 
distances ( 11 . 14}) and ( 11.151) . In momentum space, the smearing (12.131) retains only small 
wavenumbers Q, by which we mean Q whose components are all < l/£ smea r- If we make an 
approximation to (I2.12p that changes the integrand for large Q but not for small Q (where 
Q is conjugate to the point x where we measure the charge density), then the smeared 



4 Readers familiar with the (r,a) formalism may know these retarded Green functions as G ar (<i,<), 
G aar (ii, £2, t), etc. See, for example, the review of notation in Ref. [16j . 

5 Hydrodynamic quantities like viscosity and charge diffusion constants can be studied using two-point 



12, 



141, Q because one may measure them by studying the low-wavenumber response of the 



correlators 

system to a low-wavenumber source. The reason one has to go to 3-point functions here (and in Hofman 
and Maldacena 



15]) is that we are specifically interested in a high- momentum source in order to study 



"jets. 
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response ( 12 . 13[) containing all the information we are interested in will not change. In the 
rest of this paper, we will not again explicitly reference the smeared response f !2. 13j) . but we 
will feel free to make approximations that are only valid when the components of Q have 
magnitudes small compared to E and T. Note that this only applies to the Q conjugate 
to the measurement point x; no such approximation would be acceptable for the source 
momenta Qi and Q 2 in (I2.12p . which are both large. 

In particular, using the explicit source (II. 2p in ( 12.1 2p gives 

Ml [ e a e p G£ +3)aPtM (Qi, Q2; Q) A* L (Qi - k) A* L (Q 2 + k) e~ iQ ^e~^ 

JQaQ.i 



-Q1-Q2 
(2.14) 



where we have ignored terms involving A* L (Qi — k) A* L (Q 2 — k) and A* L (Qi + k) A* L (Q 2 + k) 
because these contribute only to to very large momenta Q = —Qx — Q 2 — ±2/c. 

Before moving on to the gravity side of the calculation of retarded correlators, it will 
be useful to review the fact that retarded real-time correlators are related to time-ordered 
imaginary-time correlators by analytic continuation in frequency. For two-point correlators, 
we are used to seeing this in the form 6 

G R » = G( W ±«), (2.15) 

where G is the analytic continuation of the imaginary-time Green function Ge to real- 
time frequencies; 7 the upper and lower signs are for the retarded (R) and advanced (A) 
Green function respectively; and we have suppressed showing the spatial momentum q. To 
understand the generalization to n-point functions, it is useful to write the 2-point function 
in terms of two momenta, trivially related by momentum conservation: 

G(Qi; Q) = J d A Xl d 4 x G(x x \ x) e ~ iQl ' xl e~ iQ ' x = G(Q) (2n) 4 5(Q 1 + Q). (2.16) 

Since Qi = —Q, the prescription ( 12 . 1 5[) is in this language 

G Rt x(u 1 ;u) = G(u 1 ^ie;u±ie). (2.17) 

That is, the frequency associated with the response has a ±ie prescription and that asso- 
ciated with the source has the opposite. This generalizes to the higher-point functions, so 
that [l8| 8 

Gn(u>i , oj 2 ] u) = G(oj\ — ie, u 2 — ie; u + 2ie). (2-18) 



6 The imaginary-time n-point Green function is defined here as (— )™ 1 times the imaginary-time time- 
ordered correlator of fields. So, for instance, Ge(?) = — 1/(<Z 2 + rn 2 ) for a free massless scalar. 

7 u+ie (uj—ie) indicates that one continues from positive (negative) imaginary frequencies. 

8 One may check that the ie prescriptions (|2.18[) enforce vanishing of the retarded correlator GrR,(ti,ti.;t) 
unless t is the largest of the three times, just as in the equivalent but more explicit formula (|2.1ip . (i) 
Fourier transform (|2.18|) back to t\, £2, and t, (ii) use frequency conservation 5(u>i + 102 + w) for the 
oj integration to rewrite e - lLJ i t i e - l "2t2 e -jwt ag e -iu>i(ti-i) e -iw 2 (t2-t) ^ an( ^ close the uj\ and/or LU2 

integration contours in the lower half plane if £1 > t and/or ti > t [realizing that the prescription (|2.18|) 
makes Gr(wi, W2; ui) analytic in the lower half planes of wi and/or ui2 (corresponding to the upper-half 
plane of u> = —u>i — ^2)]- 
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boundary u=0 

(a) (b) 

FIG. 4: Witten diagram for (a) 3-point boundary correlator in imaginary-time AdSs-Schwarzschild 
and (b) retarded 3-point boundary correlator Gr(xi,X2',x) in real-time AdSs-Schwarzschild. The 
darker shaded region shows the region of bulk vertex position that gives non- vanishing contribution 
to the retarded correlator, which is the intersection of the causal future of x\, the causal future 
of X2, and the causal past of x. We have taken artistic license when drawing the boundary of the 



Penrose diagram with all four sides straight 191 ] . 



C. Gravity Dual: 3-point functions 

1. Basics 

The AdS/CFT correspondence translates the problem of computing the 3-point Green 
function of R charge currents in strongly-coupled gauge theory to the problem of computing 
3-point boundary correlators of classical gauge fields living in AdSs-Schwarzschild space. If 
we were interested in imaginary-time correlators, the boundary correlator would be given by 
the Witten diagram of fig. Hk, with the circle representing the boundary of imaginary-time 
AdS 5 -Schwarzschild space. 

To start with a simpler example, if we were studying a 3-point boundary correlator of 
5-dimensional scalar fields in the gravity theory with 3-point vertex A, then the Witten 
diagram would give 

G E (Qi,Q2,Q 3 ) = A / du^g E {Q 1 ,u)g E {Q 2 ,u)g E {Q 3 ,u), (2.19) 

Ju B 

where the Qi are 4-momenta in the boundary theory, u is the coordinate for the 5th di- 
mension, and Ge{Q,u) is the bulk-to-boundary propagator which solves the 5-dimensional 
imaginary-time classical equation of motion for the scalar field, appropriately normalized on 
the boundary, ub — > is the usual boundary regulator. In our calculation, there will turn 
out to be no divergences as u — > 0, and so we will set ub = for simplicity. 

Imaginary-time AdSs-Schwarzschild spacetime is smooth, and we can integrate over the 
entire space-time without worrying about horizons or singularities. In real time, one in 
principle has to worry about such issues, but for the retarded 3-point propagator one can 
find the correct prescription simply by the analytic continuation (I2.18P of the imaginary-time 
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result: 



Gr(Qi, Q 2 ; Q) = A / duy/^g G{(jJ X - ie, Q u u) Q(co 2 - ie, Q 2 , u) Q(uj + 2ie, Q, u) 
Jo 

= xf duV^g A (Q 1 ,u)g A (Q 2 ,u)g R (Q,u), (2.20) 

Jo 

where Q R and Q A solve the linearized 5- dimensional real-time equations of motion with 
retarded or advanced boundary conditions respectively. The result (I2.20p was found in Ref. 
(20| . where more discussion of both retarded and other 3-point correlators may be found. 



(See also the related discussion in Ref. |21|.) 

As a matter of convention, note that our 4-momenta Qi, Q 2 and Q are momenta in 
the gauge theory, and therefore they are the momenta conjugate to the boundary points 
in the bulk-to-boundary propagators. (The momenta conjugate to the 4-position of the 
bulk point are correspondingly —Qi, —Q 2 , and —Q.) As a result, our convention is that a 
retarded bulk-to-boundary propagator refers to the case where information flows from the 
bulk to the boundary, and so corresponds to a solution where waves flow out of the horizon. 
Similarly, our advanced bulk-to-boundary propagator is the solution where waves flow into 
the horizon. 9 The two are related by 

Ga(q,u) = [g R (q,u)}* = G R (-q,u). (2.21) 

Readers who find it awkward or confusing to think in terms of bulk-to-boundary flow rather 
than boundary-to-bulk flow may, if desired, rewrite our equations in terms of boundary-to- 
bulk propagators <S(p,u) defined by 

®r(p,u) = G A (-P,u), & A {p,u) = g R (-p,u), (2.22) 

where p is the 4-momentum conjugate to the bulk position. 

The real-time Witten diagram associated with (I2.2(jp corresponds to Figure Sb, and the 
bulk vertex lives exclusively in the right-hand quadrant of the Penrose diagram. This can 
be understood from the causality properties of the retarded and advanced bulk-to-boundary 
propagators in the analytic continuation (I2.20p of the imaginary-time result (12.191) . Q A only 
has support when the bulk point is in the causal future of the boundary point, and Q R 
only when the bulk point is in the causal past of the boundary point. Taking the boundary 
points to all be on the right-hand boundary, the combination of these causality constraints 
requires the bulk vertex to be in the right-hand quadrant. This argument is similar to 
the discussion by Caron-Huot and Saremi [22| of the "causal diamond" in their analysis 



9 Let us relate this to the notation of Son and Starinets 12} , Let q be the 4-momentum conjugate 
to the boundary position and p = —q the 4-momentum conjugate to the bulk position. Then, in the 
scalar case, our <7r(q, u) corresponds to their f£(u). Note that /* = fl q = f q . In position space, our 
retarded bulk-to-boundary propagator from a bulk point (y,u) to a boundary point x is Qb,(v,u;x) = 
Iq fq( u ) e %q ' iyX ~ v ^ ) — J fp(u) e lp '( y ~ x \ The last is the same as the advanced boundary-to-bulk propagator 
&A(x;y,u) defined by (I2.22p . If one wished to describe a bulk wave falling from the boundary into the 
horizon, it would be given by <p(y, u) — J d 4 x ©R,(a;; y, u) <fi(x, 0) = J d 4 x GA(y, u; x) <j>(x, 0) or equivalently 
<t>(p, u) = R (p, u) (f>( P , 0) - g A (-p, u) <t>(p, 0) = g A (q, u) <t>*{q, 0). 
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of long-time hydro dynamic tails from one-loop gravity corrections to the 2-point retarded 
correlator. 

It is the simple form (12.201) of the retarded correlator that led us to choose to set up 
the problem of studying jet evolution in terms of retarded correlators, rather than the type 
of correlators (0^ q O(x) O q ) considered by Hofman and Maldacena, which would be more 
complicated to analyze at finite temperature. 10 



2. The vector-vector-vector vertex 

The generalization of (I2.2(jp from scalar interaction to vector interaction is that we need 
to use the vector-vector-vector vertex in the super-gravity (SUGRA) theory: 

(Gr) { £ = /^«V^vertex^ c )^(^ A a (Q 1)W ),^(Q 2 ,«),^ M (g )W )), (2.23) 
J o 

where the advanced and retarded bulk-to-boundary propagators Qi a have one 4-dimensional 
vector index (a) associated with the boundary point and one 5-dimensional vector index (I) 
associated with the bulk point and are normalized on the boundary so that 

Gp a (Q,0)=VPa and G 5a {Q,0) = 0. (2.24) 
The vertex function is extracted from the cubic terms in the SUGRA interaction (si, Hi| 

- / d 5 x^F IJ «F?j ~^- 2 Jd 5 x [d abc e IJKLM AUdjA b K )(d L A c M ) + ■■•], (2.25) 

where we work in real time, have only shown explicitly the cubic term in the Chern-Simons 
interaction, and 

9sg = j^ and k = N 2 c -l. (2.26) 

The F 2 term produces a 3-point vertex with R charge structure f abc , whereas the Chern- 
Simons term (which reproduces the R charge current anomaly) yields d abc . Here f abc and 



10 Such correlators can be related to the retarded correlator. Since the source is localized, and we are 
only interested in measurements at times x° after the source turns off, we can rewrite the correlator 
(Oj 0(x) O q ) in Schwinger-Keldysh (closed time path) notation as G2ii(Oj, O, O q ). Each subscript in 
G211 designates whether the corresponding operator appears on (1) the first leg of the Schwinger-Keldysh 
integration contour (t—— 00 to +00) or (2) the second leg (+00 to —00), and operators are ordered 
accordingly. Specifically, G^ = (A(J~BC)) , where T represents ordinary time ordering. (Since a; is 
our largest time, G211 = G221 = G2 r i here, where r is the average of using 1 and 2.) 3-point Schwinger- 
Keldysh propagators G aia2a3 can all be (non- locally) related to retarded Green functions G aar , G ara , and 
G raa (the distinction being which of the three operators is the latest time when defining "retarded") and 



their complex conjugates. Explicit formulas (and more explanation of the notation) are given in Ref. 16|. 



Alternatively, see Ref. 20( for a discussion of different 3-point correlators directly in terms of integrating 



over both right and left quadrants in the gravity theory. (Beware that Ref. 20( , following Ref. 23j , uses 



a different convention for the Schwinger-Keldysh contour than Ref. 16| . This introduces factors of e 
into the definition of G aia2a3 , where a = 0/2.) 
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d abc are defined in terms of SU(4) Hermitian generators T a by tr(T a T b T c ) = \{d abc + if abc ), 
normalized by tr(T a T b ) = \ 5 ab . 

Recall that the currents in our source fl X . 2 1) and measurement fl 1.7ft were chosen to lie in 
an SU(2) subgroup of SU(4). Since d abc vanishes for SU(2), the Chern-Simons term will not 
contribute in our application. 11 So the only relevant contribution to the vertex comes from 
the cubic interaction 

yafec 



[ d 5 xV^gg IM 9 JN (d I A a J - djA?)A h M A%. (2.27) 

/SG-"- J 

The factors of the AdS radius R all cancel [the explicit i? _1 above with the factors in 
^/—gg IM g JN from (12. ip ]. If desired, one could simply set R — 1 in the rest of the paper. 

To get the vertex function in (12.231) . substitute the three bulk-to-boundary propagators 
G for the three A's in (I2.27P in all possible permutations. Since the problem studied in this 
paper is invariant under translation in the transverse directions, we will restrict attention to 
5-dimensional gauge choices that respect this invariance. In particular, the bulk-to-boundary 
propagators will preserve transversality of polarization, 

e a Gi* oc ej, (2.28) 

and transverse derivatives will vanish, 

e a g aI dj ■ ■ ■ = 0. (2.29) 
Here it is convenient to define a 5-dimensional ii in terms of 4-dimensional e a by 

£j = (e , e 1} e 2 , e 3 , 0) = (0, 1, 0, 0, 0). (2.30) 
Putting everything together, the piece of (I2.23P that contributes in our problem is then 



jabc 



j dV du V=gg IM g JN {l-eid'j GtW, % Xl )] e M g^x', u- x 2 ) g^(x' 



x , u: x) 



-{x 1 ^x 2 )\, (2.31) 



where G± = e p Q pa e a = £ IJ ,'r] w GpaTf u E lJ * s the transverse piece of the bulk-to-boundary prop- 
agator. (Other than as an attempt to save space in equations, there is no significance to 
whether we write the R/A for retarded/advanced as subscripts or superscripts.) Switching 

to four-dimensional momentum space, and using the notation fdg = fdg — (df)g, 

£abc r 

^(GWX? = -is I duV^g(eig IJ ej) 



gsc R 

X 



GiiQi, u) £tf (Q 2 , u) t(Q 2 - Qi) p! rg*(Q, u) 



+ QKQuu) d 5 £tf(Q 2 , u) g 55 Gl(Q, u) 



(2.32) 



An independent reason that the Chern-Simons term does not contribute in our problem is our choice of 
linear polarization e for our source: the two A's which attach to the source points on the boundary will 
give factors of e and e* oc e contracted with the s IJKLM , giving zero. 
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3. Summary in ^4s=0 gauge 



The last formula fj2.32j) is simplest in v4s=0 gauge, where the last term vanishes. We will 
specialize to A 5 =0 gauge in the remainder of the paper. Define G^ ± by 

^(Gr)<S? = f abc Gl^. (2.33) 
Putting everything together and using / 1-3 = 2i, we then have 

~ Wsftr I GliAQu Q2; Q) A* L (Qi - k) AKQ2 + k) e^-e***-* 

JQaQ.i 



(2.34a) 



with 

<*t±p = KQ2 - Qi) P -gTjiJ du ^9 (ei9 IJ ej) sPQi(Qi,u) Gt(Q2, u) f?» (Q, u) (2.34b) 
in A 5 =0 gauge. 

As it currently stands, ( I2.34p is challenging to evaluate. Our crucial approximation in 
what follows will be to replace G R (Q, u) in (I2.34bl) by its small-Q approximation valid for 
Q < T. As discussed earlier in section III B\ such approximations are adequate for the 
long-distance physics that we wish to study. As a simple example, for a massless bulk scalar 



field, the bulk-to-boundary propagator in the small-Q limit is [12|, [25 
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QLUQM = (1 -- 2 )"^ T + 0(^,i^), (2.35) 
where u = Q°. The approximation (I2.35P has the nice property that it factorizes for Q = 

Ql*~{Q> «) ^ (1 - uJ^ T {\ - uJ^ T . (2.36) 
If this were the propagator to use for Q R (Q,u) in (|2.34bl) . then the Qi and Q2 integrals 
would factorize in (I2.34ap . greatly simplifying the calculation by allowing us to evaluate 
them independently. We will see in section [IV] that the issue of factorization is a little more 
complicated for the vector bulk-to-boundary propagator G^AQi u ) a ^ finite temperature, but 
we will be able to use a variant of this trick to factorize the calculation. 



III. THE ZERO-TEMPERATURE CALCULATION 



In this section, we warm up to the calculation by applying to the case of zero temperature 
the methods that we will later use for finite temperature. At zero temperature, the metric 
(12. ip becomes 

'1 , .0 . 1 



ds 1 



4 



-dt z + dx z ) + 
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u 2 



du 



(3.1) 



12 To have an expression valid all the way to the horizon, it is important not to expand in powers of the 
exponent as (1 — u 2 )-"*>/ 47rT ~ 1 — ln(l — u 2 ). For any fixed oj <C T, the corrections to this truncation 
are not small when u is close enough to the horizon that (oj/T) ln(l — u) 3> 1. 
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for AdSs, where u runs from zero (at the boundary) to infinity. This is related to the 
Poincare metric 

ds 2 = R 2 ^ dx ^ dx - ± dzl (3.2) 
z 1 

by u = \ z 2 . Eq. f!2.34bj) for the 3-point function becomes 

1 /*00 7 — 

Gl ± , = ^(Q2-Ql) P V p(T V ^ / —G^Qi,u)g^Q 2 ,u)G^(Q,u). (3.3) 

At zero temperature, it is notationally a little more convenient to use the coordinate z 
than u. However, we will use the coordinate u for the finite-temperature calculations, and 
so we stick to u here in order to make the two cases as easy to compare as possible. 

In imaginary time, the regular solution to the linearized 5-dimensional vector equation 
of motion d *F = is 

Ap(q,u) = Gpv{q,u)A v (q,0) (3.4a) 
in momentum space in ^5=0 gauge, where 13 

G*Mu) = vW^i(vW) (V - q ~f) + ^ (3-4b) 

and K n is the modified Bessel function of the second kind. The real-time version is then 
simply 

vW K x (vW) U„ ~ qj f)+ qj f- (3-5) 

where q 2 = r] al3 q a qp with u — > u> ± it for retarded and advanced. For time-like momenta 
(q 2 < 0), it can be useful to recast K\ in terms of Hankel functions, 



/7~T is l ^F _ 2^ / in V '~ U< 1 2 H i ] '(V '-^m 2 ) for retarded q° > 0; 

VW^l(VW) = _. _ r-— 2 „( 2 ) _ r—A . n ^ . (3.6) 



A. The crucial approximation 

Now we will make the same approximation that we suggested for the finite temperature 
limit: take the small-Q approximation of G^iQ^u). At zero temperature, however, we do 
not have temperature to define a natural scale Q < T, and so the discussion of the small-Q 
limit is a bit more complicated. The zero-temperature G atl {Q,u) given by (I3.5P simplifies 
when uQ 2 <C 1, in which limit it is simply r] afl . 

The crucial approximation: 

0*(Q,tZ)^7^, (3.7) 
where Q = —Qi — Q2 is the momentum conjugate to the measurement point x. 



13 



For the momentum-space solution in covariant conformal gauge (V 1 Aj = 0), see, for example, Ref. [2 
specialized to the case d=A and m—0. This is the gauge transformation Aj —> Aj — djA of (|3.4p with 
A(q, u) =l[2uq 2 K 2 {y/titf) - 1] (q u /q 2 ) A u {q, 0). 
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We will discuss the validity of this approximation in a moment, but first let's see what it 
gives. Making the zero-temperature version (13. 7p of the crucial approximation in (13. 3 p and 
using the relation (" 12 . 2 1 [) between advanced and retarded propagators, we may approximate 
fICTaj) and (E3D as 

{f(x)) ~ -rT / —A(x,uyid p A(x,u), (3.8a) 
9sg Jo u 



where 



A(x, u) = [ gf(q, u) A L (q - k) e iq - x . (3.8b) 

J a 



(Note: q ■ x is a flat-space dot product rf v q^Xy here.) The combination —A*id p A has the 
same form as the expression for the "current" associated with a charged bosonic field in field 
theory. 

The definition (13.8bj) of A is simply (up to normalization factors) a convolution of the 
retarded bulk-to-boundary propagator and the positive-energy piece of the classical source 
(11. 2p in the 4-dimensional gauge theory: 

ri,^A(x, u) cx J d A y g* a (y, u; x) A« l(+) (y) (3.9) 

where 

Al [+) {y) = ^Na^ K L {yl (3.10) 

That is, A(x, u) is a retarded solution to the linearized 5-dimensional vector equation of 
motion that is proportional to the source A c u + ) on the boundary. As time x° progresses, the 
bulk excitation represented by A will fall away from the boundary, into the fifth dimension. 14 
This evolution of A is roughly the type of problem studied in the context of jet evolution 
by Hatta, Iancu, and Mueller [5(. 15 They made qualitative interpretations concerning the 
initial conditions on A and the result of its evolution. Here, we have seen the question of 
.A's time evolution arise step by step from a problem posed completely in the 4-dimensional 
field theory. Consequently, we also have a quantitative way to interpret the solution for A: 
use it in (I3.8ap [or the appropriate generalization to finite temperature coming in section 
UVj to find the current density response. 

When does one expect the approximation (13. 7p to be valid in the zero temperature case? 
It is correct in the limit FPg^Q^Qu = AuQ 2 < 1, where Q 2 means rf v Q p Q v = AQ + Q_. 



4 Our nomenclature gets a bit convoluted here. Earlier we said that our convention was that Q corresponds 
to waves propagating from the bulk to the boundary. So why do we say here that A(x, u) falls into the 
horizon as a; increases? Note that x is the boundary point for G^iy, u; x), whereas (y, u) is the bulk point. 
x° — y° — > oo is equivalent to y° — x° — > — oo: the bulk point becomes further and further back in time 
compared to the boundary point and so closer to the horizon. It might have been clearer in this respect 
to use notation like A(u; x) or A(u — > x) instead of A{x,u) [and perhaps similarly for G K {q, u)], but we 
decided that would be too cumbersome. 

5 A very minor difference is that for most of their paper they chose to focus on the evolution of the Aq 
and A3 components of the 5-dimensional field, whereas we have found it convenient to choose a transverse 
source, and so our A tracks the evolution of the transverse components. 
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Writing Q± ~ 1/Ax where Ax ± are the desired resolutions of x + and x in the response, 
this suggest that the approximation is valid when 

Ai + Ai"»m. (3.11) 

The question then boils down to knowing the natural scale for u in this calculation. 

Later on we'll see that one of the natural scales for u that arises in the analysis, identified 
by Hatta et al. [5|, is 

u~^. (3.12) 
However, we will find that a more important and parametrically larger scale is 

(3.13) 

The approximation (13. 7p is then valid provided we only apply the result to resolve questions 
on scales 

(x + ) 2 

Ax+Ax->^— ^. (3.14) 
EL 

Since EL 3> 1, there is no problem simultaneously resolving Ax + and Ax~ to scales small 
compared to x + . 

If we want to resolve Ax~ all the way down to Ax _ ~L (as we shall implicitly do later) and 
simultaneously resolve Ax + to better than x + itself, then (13.141) requires x + <C EL 2 . These 
constraints are special to the zero-temperature problem — at finite temperature, choosing a 
resolution distance large compared to 1/T will be all that we will need. 

Thanks to the approximation (13. 7p . we now have just a single 4-momentum integral 
( 13.8bj) to evaluate or approximate, followed by a one-dimensional u integration (I3.8aj) . In 
Appendix |XJ we verify that the approximation ( I3.8a|) obeys current conservation outside of 
the space-time region of the external source. 



B. The high-energy approximation 

Our next task is to evaluate A. Using q 2 = ^q+q~ and rewriting g_ as E + Ag_, the 
formula (I3.8bj) for A becomes 



A(x, u) = e iE *~ J 2 Gl{±u{E + Ag_)g + ) Ag>( g+ , Ag_) e ^ + e ^^^ ( 3 .15) 

where the transverse bulk-to-boundary propagator extracted from (13. 5p depends only on the 
combination uq 2 : 

g ± (uq 2 ) = Q±(q,u) = VWifiiVS?). (3.16) 

A L (q + , q_) is the two-dimensional Fourier transform (II. 18ft of the source envelope, which 
for a Gaussian envelope (11.61) would give 

Ai 2) (g + , Ag_) = 2ixL 2 e-^ 2 e-^- L ?. (3.17) 

The smooth envelope function restricts support for the integral (I3.15P to Ag_ < 1/L and 
therefore Aq_ -C E. We may therefore approximate 

Gl(Au{E + Ag_)g+) ~ gf(4uEq + ) (3.18a) 
19 



x<> x+ 




x~ 

FIG. 5: The space-time development of a high-energy excitation at zero temperature. 



in (EI5]). 

With this approximation, the Ag_ integration gives 



A(x,u) ~ e iEx ~ J ^ gf(4uEq + ) A^(g + ; x~) e iq + x+ , (3.18b) 



where 



kf{q + -x-) = J dx + A L (x)e- tq+x+ (3.19) 
is insignificant unless q + < 1/L and x~ < L. For the Gaussian envelope ( II. 6p . 

Ai 2) (g + ;x") = 2 7 r 1 / 2 Le-^ i >V^) 2 / 4L2 . (3.20) 

As depicted in fig. [5j the response (13.18bj) is localized in x~ . Up to the matters of 
resolution discussed at the end of section IIII Al the high energy excitation produced by the 
source simply propagates along the lightcone to the right of the 1+1 dimensional source 
region. 

Now return to eq. ( I3.8aj) for the current response: 

K[2 poo i- 

(f(x)) ~ - v w —A(x,uYid p A(x,u), (3.21) 

#SG Jo u 

The derivative d p A of (13. 18b[) will be dominated by the term where the derivative hits the 
highly-oscillating factor of e lEx : 

(f( x ))^2k» J -^ —\A(x,u)\ 2 . (3.22) 

9sG JO U 

In this approximation, the statement d^j^ = of current conservation (outside of the 
source region) is d + (j + ) = 0, which means that (for zero temperature) 



00 dfi 

— \A(x,u)\ 2 (3.23) 
u 



should be independent of x + outside of the source region. 
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x<> x+ 




x- 



FIG. 6: Picture of the source and response after the (flawed) approximation ()3.25j) or ()3.26p . 



C. An approximation that doesn't quite work 

We now discuss a further approximation, which will be flawed. But it will be close to 
what we eventually need, and how it fails will be instructive. We will fix it up afterward. 
In the case of a Gaussian envelope (I3.20p . the evolution f!3.18bj) of A is 

A(x, u) ~ 2n l ' 2 L e iEx ~ e'^ 2 '^ [ ^± gf(AuEq + ) e~ (q+L)2 e iq+x+ . (3.24) 

J 27T 

The integral determines the x + dependence of the response. There are three factors in the 
integrand, associated with three different scales for q + : the corresponding wavelengths l/q+ 
are 4uE, L, and x + respectively. To study the response far away from the source, where 
L C x + , one might hope one could treat L as arbitrarily small in the integrand of ( 13.24}) . 
replacing it by a convergence factor: 



where e is infinitesimal. Mathematically, this approximation corresponds to replacing the 
Gaussian source envelope 

A L ( X ) = e -\{* + /L? e -\(x-/Lf _^ 2^ 2 Le-^ x+ / L ^5(x-) (3.25b) 

and so corresponds to replacing fig. El by fig. El The generalization to generic source envelopes 
would be 

A(x, u) -> e iEx ~ Af (0; x~) [ — Gf(AuEq + ) e~ eq2 +e iq+x+ . (3.26) 

J 27T 

Changing integration variables from q + to 

K = AuEq + , (3.27) 

we can rewrite (13. 26j) as 

Mx,u)^e^-^fl,(^-), (3.28) 



AuE \AuE 
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where 

I(s) = J d ^Ql(K)e-^e^. (3.29) 

Note that the only parameter in the definition of I(s) is its argument s, and so a natural 
scale of the problem is s ~ 1. From (I3.28p . this corresponds to 

u ~ ^, (3.30) 
which was the first of the two scales for u previewed back in (13. 12ft . 



An apparent paradox 

We can now see the origin of a problem. The response (13.221) after all our approximations 
becomes 

oL/i Kf2 poo 

^{x))^-— J -^\kf(V-x-)\ 2 dss\I(s)\ 2 , (3.31) 
\ x ) 9sg Jo 

where 

s - 4 < 3 - 32 > 

is oo on the boundary u—0. This result appears to depend on x + as l/(x + ) 2 , which is 
inconsistent with current conservation, as discussed for (13.231) . The loophole to this paradox 
is that the integral in (I3.3ip turns out to be s— >0 divergent (corresponding to contributions 
far away from the boundary). Specifically, explicit evaluation of the integral J(s) defined by 
(EEHgp yields (see Appendix |BJ) 



I(a) = ~e i /'B( 8 ), (3.33) 
where 9(s) is the step function, giving 



s 2 



dss\I{s)\ 2 oc I ^. (3.34) 



As we will now see, this divergence arises from ignoring the width L of the source in x + in 
order to make the last approximation (I3.26p . 



D. Fixing the last approximation 

We will now step away from the last approximation (I3.26P but will still treat L as relatively 
"small" in a sense to be made precise in a moment. We know from the previous analysis 
that what will be important are small values of s = x + /4uE. Since the full integral in eq. 
(I3.18bl) for A is more complicated than I(s), we will make our lives easier by approximating 
it in the small s (large u) limit instead of attempting to evaluate it exactly for general s. As 
we'll see, the appropriate approximation in this limit is the method of steepest descent. 
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1. Steepest descent evaluation of A 



For the sake of concreteness, we focus for the moment on the case of the Gaussian envelope 
function, which gives f)3.24p for A, proportional to 

X(x + ,u) = f ^±gf(4uEq+)e-^ 2 e iq+x+ . (3.35) 

J 27T 

The first thing to note is that the contribution to this integral from space-like momenta 
(g + >0) is small in the limit u — » oo. That's because, in that case, Q^{AuEq + ) falls rapidly 
with q + . The dominant large u behavior therefore comes from the time-like (g + <0) region 
of integration, where Q R is oscillatory. Writing q + = —k — ie (retarded prescription), the 
integral then gives 

roc a-i 

l(x + ,u)~ / —g ± {e^AuEk)e' {kL)2 e- ikx+ . (3.36) 
Jo 2n 

Now use the large argument approximation to eq. (13.161) for Q± (which we will verify later 
is appropriate): 

g ± {£) ~ v^^V 2 ^, (3.37) 

giving 

l(x + , u) ~ e-™'\k (\QuE f A / — kWe- 3 ™, (3.38) 
v Jo 2?r 

where 

S(k) = -iVl6uEk + ikx + + (kL) 2 . (3.39) 

We now perform steepest descent by finding the zero k=k+ of dS/dk. Formally expanding 
the solution in powers of L and keeping only terms through L 2 , this extremum is at 



AuE 



x+) 2 



1 + l T^ L 



(3.40) 



The condition for treating L as small in this way is that the magnitude of the L 2 term be 
small compared to that of the L° term, which is the condition 

(x+) 3 (x+) 2 x + , s 

u<< -eW--eTT- ^ 

The value of S corresponding to (13.401) is 

x 4uE UuEL) 2 



a; T (x n 



giving 



e s ^*^ ~ e i/s exp 



(AuELf 
(x+Y 



(3.43) 



The e l l s factor is just the oscillatory factor in the earlier result ( I3.33|) for lis). The size 
|«5(Av)| of the exponent will be large, justifying the steepest descent approximation, when s 
is small. By considering the first correction to the exponent in powers of L, we have now 
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found in ( I3.43|) that there is a decreasing exponential that kicks in and cuts off the large u 
behavior. This large-w suppression factor starts to kick in when 

(x + ) 2 

u ~ VLL, (3.44) 

which is the second scale for u previewed in f 1 3 . X 3 j) . At this u, the small- L approximation 
(I3.4ip is valid (far away from the source, x + ^> L), and so our treatment of L 2 effects as a 
perturbation to the evaluation of iS(fc*) is justified. 

The reason we could not just stick with the leading order term —i/s in the exponent 
S(k+) is because it was pure imaginary: we had to go to first order in L 2 to find the leading 
contribution to the real part of the exponent. In contrast, a leading-order evaluation of the 
algebraic prefactor of the exponential is good enough, and so we can ignore the effects of 
L on that prefactor. The upshot is that (I3.43P tells us to modify our previous analysis of 
section IIII CI by simply replacing 



I(s) — > I(s) exp 



{AuELf 



3T 1 



(3.45) 



Alternatively, one may obtain the same result at small s by explicitly finishing the steepest 
descent analysis by expanding S{k) to quadratic order about S{k*) and doing the Gaussian 
integral to get the prefactor. 

From (I3.24p . (13 .33 p . and (I3.45p . the final expression for A is 



A(x, u) ~ -i2-K 1 ' 2 e ^- e -(^) 2 /4^ e mE/ x + exp 

[x + y 



(AuEL) 
(x+Y 



6(x + ). (3.46a) 



The generalization of the above analysis to generic source envelopes simply replaces the 

(2) 

large-w suppression factor above by the envelope function A L (q + ;x~) evaluated at the 
saddle point q + = —k*: 

A(x, u)^-i^ e—/* + Ai 2) (-^; *") °(* + )- (3.46b) 
2. Final result for current response 



Inserting f!3.46bp for A into the expression (I3.22p for (j M (a;)), and changing integration 
variable from u to q + = — AuE/(x + ) 2 , gives the final zero-temperature result 

\[2 fO 

(f( x ))^2k^9(x + ) / dq + \q + \ \A^(q + ; x~)\ 2 (3.47a) 

for |a; + | ^> L, with g$Q given by (I2.26p . This result is independent of x + after the source turns 
off, as required by the discussion of current conservation surrounding (I3.23P [which we remind 
the reader is subject to the resolution requirement (I3.14p for making our approximations]. 
For the case of the Gaussian envelope (13.201) . the result (I3.47al) is 

\[2 

(f(x)) ~ 2<kW e ~ (x ~ )2/2L2 6(x + ). (3.47b) 

S'SG 
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The total charge density per unit transverse area at late times can be taken from (I3.47aj) 



as 



Q 



dx 3 (j°(x)) 

AT 2 f f° 
2E ^ / dx~ / ^+|g + ||A[ 2) (g + ;x-)| 



which may be rewritten as 



Q ~ 8nE 



Ma f 2 dq + dq. 



e(-q + )\q + \\A^(q + ,q- 



(27T) 



(3.48) 



(3.49) 



[from where one may now see where the form of the normalizing denominator in f 1 1 . 1 7 j) comes 
from]. In the case of the Gaussian envelope f)3.17p . 



Q = (2nf 2 EL ^ 



(3.50) 



As a simple check of our calculation, in Appendix IA2I we make a completely independent 
computation of the total charge produced by the source by using Ward identities. The result 
is the same as (13.491) . 



E. Schrodinger interpretation d la Hatta, Iancu, and Mueller 



In 4-momentum space, the zero-temperature, linearized equation of motion d*F — for 
the transverse components of the 5-dimensional vector field is 



u 



0. 



(3.51) 



The solution, appropriately normalized, is just the transverse bulk-to-boundary propagator 
(13.16!) . If we approximate q 2 ~ 4Eq + and Fourier transform q + to x + , the above equation 
can be rewritten as 

/ i r \ 

(3.52) 



Our solution (13.461) for A(x,u) (approximately) solves (I3.52p . 16 Focusing on the Gaussian 
source case for concreteness, the x + and u dependence of A is 



A oc MM exp 

( X ) 



{AuELf 



lx n 



(3.53) 



Hatta, Iancu, and Mueller [5j noted that the equation (I3.52p (as well as its finite- 
temperature version in the case of small u) can be recast in the form of a Schrodinger-like 
equation by changing variables from u to z = 2^/ti and redefining 



A±(x + ,ll) = y/z~(j)(x + ,z). 



(3.54) 



16 By direct substitution of (|3.53j) in (|3.52[) . one may double-check that it is an approximate solution 
everywhere in the region u < (x + ) 2 /EL, which is the region where A is non-negligible. 
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1 ■■ 




x+/E ~ {x+Y/EL 




FIG. 7: Qualitative pictures of the real or imaginary parts of (a) A versus u and (b) <f> = A/yfz 
versus z = 



The equation of motion (I3.52p then becomes 



2id, 



2E + 8Ez 2 



(3.55) 



which is a Schrodinger equation with potential energy oc z 2 , provided x + is interpreted as 
time. The conserved "probability" of this Schrodinger equation is 



dz 



\ I f \A\*. 



(3.56) 



From our discussion surrounding (I3.23p . we see that, for zero temperature, conservation of 
probability in this Schrodinger problem is equivalent to conservation of charge in the original 
field theory problem. 

Qualitative sketches of our result (13.531) for A and are given in fig. [71 Both (i) the 
wavelength z r- 



\J x + / E of the first oscillation and (ii) the location z ~ x + /yEL of the 
bulk of the probability grow with time (x + ). The substantive difference with the study of 
Hatta et al. jH[ is that they studied non- localized solutions, 17 



A oc 



u 



{x+y 



i&uE/x A 



(3.57) 



which do not decay at large u and have infinite normalization J dz \<fi\ 2 . These solutions cor- 
respond to taking L — > in ( I3.53P and so, in our context, making the failed approximation 
( I3.26p . Correspondingly, the only z scale that Hatta et al. identify at zero temperature is 



17 There are some other differences. The x + in our discussion plays the role of It in theirs. Similarly, our 
A + plays the role of their Aq in what follows. They focus more on the longitudinal mode A + than the 
transverse mode (|3.57l) . The two are qualitatively similar, as can be seen from our Q3.5p . except that Hatta 
et al. choose to work with a = d a A + instead of A + directly. They choose boundary conditions so that a oc 
(x + )~ 1 e l4:UE / x+ [see their (3.13)], whereas the boundary conditions determined by our type of field theory 
problem would give a different solution (in the L~^0 limit) to the same second-order equation [their (3.11), 
dropping the K 2 term, with ip and a related by their (2.8)]: a = d a A oc (x + )~ 2 (l + iAuE/x + )e lAuE l x+ 
from our f|3.57|) . 
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On 



FIG. 8: Qualitative picture of the two u scales of fig. [TJversus x + . The top of the plot represents u 
very close to the boundary (u ~ L/E); lower points are deeper in the fifth dimension. The dotted 
line represents u ~ (x + ) 3 / EL 2 , above which L can be treated as perturbatively small in a steepest 
descent analysis. 



the scale z ~ x + /E of ( 13. 12ft . associated with the first oscillation, which they call the "dif- 
fusion" distance in z. However, the bulk of the probability density is instead characterized 
by the larger scale z ~ x + js/ELoi (13T31 . which grows faster with time. 

When we go to finite temperature, we will correspondingly find two time scales for how 
long it takes A to fall into the event horizon. The time scale for the first oscillation of A 
to fall into the event horizon will turn out to be the scale ( 

2^3 J max ~ 

E^ 3 /T^ 3 of ( 11.14p . as 

found by Hatta et al. But the time scale for the bulk of the "probability" density to fall into 
the horizon will turn out to be the shorter time scale ( x 3 ) dominant ~ (EL) l / A /T of (11 . 15)) . 

To facilitate later comparison with the finite temperature case, we repackage in fig. [8] the 
qualitative information from fig. [71 The u scales marked by the dashed lines in fig. [7J are 
represented on a plot of u versus x + by the dashed curves in fig. [HJ The region parametrically 
below the top dashed curve is where A has many oscillations and steepest descent is a useful 
method of approximation. The (magenta) region parametrically below the lower dashed 
curve is where the amplitude of A is exponentially suppressed. 



IV. THE FINITE-TEMPERATURE CALCULATION 

We now apply to the finite-temperature case the methods just presented for the zero- 
temperature case. In this section, we will use units where 2irT = 1. Powers of 2ttT may be 
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restored by dimensional analysis, replacing 

a? -> (2nT)x tl , 

<f -> (27rT)-V, 
L (2ttT)L, 

£ (2ttT)- 1 E 

throughout this section. 

The AdSs-Schwarzschild metric is 



ds 2 = — 



~(-f dt 2 + dx 2 ) + — du 2 

U U 2 J 



and the basic equations (12.341) for our problem become 

~ 2iMW I GliAQu Q2; Q) A-UQt - k) A* L (Q 2 + k) e -^- e -^ 
JQ1Q2 



(4.1) 
(4.2) 
(4.3) 
(4.4) 



(4.5) 



Q= 



and 



with 



1 r 1 du 



P 2fl , SG J» 



T po Gt(Qi,u)gt(Q.2,u)g«(Q,u) 



R 



u 



T 1 



•pa 



R 2 gP a 
Au 



pa 



(4.6a) 



(4.6b) 



(4.6c) 



A. The crucial approximation 

1. The small Q limit 

The first thing we need is the low momentum approximation for the bulk-to-boundary 
propagator (?? u (Q,u) associated with the measurement point x in (14. 6p . The issue of res- 
olution scale will be much more straightforward than in the zero-temperature case: Here 
we will simply restrict attention to distance and time scales large compared to 1/T, and so 
small Q will mean that all components of Q are small compared to T. 

Because our problem is transverse-translational invariant, Q± = and we focus on u> = 
Q° and k = Q 3 . In Appendix [UJ we show that the small-Q form of the bulk-to-boundary 
propagator in A 5 =0 gauge is, to leading order in the size of u and k 2 , given by 18 

G*(u, k) n^au ~ - — — rjjfa'ao + ka 3 ) - — ^(1 - u) l -^' 2 {ka Q + ua 3 ), (4.7a) 

g*(cu, k) tTou * -—^(^0 + ka 3 ) + —^—r 2 {\ - u)- w / 2 (ka + ua 3 ), (4.7b) 
and (irrelevant here except for comparison) 



We are interested in leading order in w and k 2 , and not simply leading order in u> and k, because we are 
interested in studying diffusion, for which u ~ k 2 . 
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g%(u, k) if v a v ~ (1 - u)-^ 2 a ± , (4.7c) 

over the entire range of u. The arbitrary 4-vector a M above represents the source on the 
boundary and is introduced here as a notational device for displaying the individual com- 
ponents of Q K . Note that the above expressions have the property that 

g*(u,k)rrQ» = Q*, (4-8) 

which is a general property of the bulk-to-boundary propagator in ^5=0 gauge, arising from 
gauge invariance [see (1A17|) in appendix [A] . The small-Q form (j4.7cj) of is irrelevant 
simply because the transverse-translation invariance of our source implies that there will be 
no expectation of j 1 - (which involves transverse derivatives of the fields) in the response. 

Eqs. ( 14. 7p are a bit cumbersome to deal with, and it will greatly help to first make some 
additional approximations. We will see later that the terms with factors of (1 — u)~ w l 2 will 
never be important for small 1 — u in our evaluation of (j^x)). Because of this, we can 
approximate (1 — u)~ tu1 / 2 ~ 1 in these terms, to leading order in the size of to. Then (14. 7p 
simplifies to 

Go^u, k) if v a v ~ r^( ia o + ka 3 ) - 7^(1 - u)(ka + ua 3 ), (4.9a) 

gl(u,k) V ^a u ~a 3 . (4.9b) 

Of the terms remaining, we will later see that the one which dominates the calculation of the 
charge deposition Q(x) of (II. lip is the fist term in (I4.9ap . provided we only wish to resolve 
Q(x) on scales large compared to the source size L as in (ll.lip . For the sake of simplifying 
the presentation, we will ignore the other terms for now and replace (14. 7p by 

g$ (u>, k) rf v a v - (za + ka 3 ), (4.10a) 

gf IM (u,k)r l ^a u ^0. (4.10b) 

Once we see how the calculation works out, we will return in sec. IIVGI to see why the 
(1 — u)~ lljj l 2 factors in (14.71) and the other terms in (14. 9p are unimportant. 

Note that the relationship g^ 3 = ikg^ in the approximation (14.101) implies via (14 .6p that 

(fix^^-diifix)), (4.11) 

which is the standard relationship between current and charge densities in a diffusive process 
(in units where the diffusion constant is 1). 



2. Factorizing the calculation 

When the formulas for g^JQ,u) taken from (I4.10p are used in (I4.6p . the Q\ and Q2 
integrals do not factorize like they did in our zero-temperature calculation. That's because 
of the ioj iibj — k 2 ) factor. We can get rid of the iou — k 2 denominator by studying the charge 
deposition function 

QQ(x) = (dt- V 2 )<j°(x)> (4.12) 
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of (11. 9p instead of directly calculating the current response (j^(x)). It will be even more 
convenient to first calculate the time integral 

S e (x) = Q I dt'G(t',x) (4.13) 

J — oo 

of the charge deposition function, which is related to the charge response by 

(d t -V 2 )(f(x))=d t i:e(x). (4.14) 

In Fourier space, 

</H-^Se. (4.15) 

The factor of iui/(iu — k 2 ) above will cancel the similar factor from (14.101) so that (14.61) leads 
to an expression for Eg where the Qi and Q2 integrations factorize. Specifically, combining 
gSD, (ED}, and {USD, 



£ (x) ~ ^ [ ^-A(x,u)*id t A(x,u), (4.16a) 

gL Jo u f 



N\ f 1 du 

gi G Jo uf 

where 

A(x, u)= I Gf(q, u) A L (q - k) e iq - x . (4.16b) 



This is very similar in form to the zero-temperature expression ( 13. 8 p for (j°(x)) except for 
the important distinction that (I4.16ap gives S e (x) instead of the charge density. 

The divergence at the horizon of the factor 1/f in the integrand of ( 14. 16al) will turn out 
to be crucial for getting a physically sensible result for J2q(x), and we will later see that 
the terms of Q^(Q,u) in ( 14. 9 p that we dropped in ( I4.10p are ignorable because they do not 
generate a similar divergent factor as u—>l. 



B. What will £©(a;) look like? 

Before we discuss the calculation of £e(x), it will be helpful to have in advance a qual- 
itative picture of what the result should look like. In fig. ^jp, we gave a pictorial rep- 
resentation of the final result (II. lip that we will find for the charge deposition function 
QQ(x) = <9 t £e(x) when resolved on scales large compared to L. The picture is that the 
excitation initially moves ballistically at the speed of light, and no charge is deposited until 
the jet reaches its stopping distance, which we will find is stretched out between the scales 
(EL) 1 / 4 and £ 1/3 . But now consider the case where we choose L>T (but still small com- 
pared to E 1 ' 3 ), and consider what happens if, unlike fig. [3h>, we resolve the charge density 
and <2© = (d t — V 2 ) (j°(x)) down to the scale L itself. At early times, before the earliest 
stopping time scale (EL) 1 / 4 , the charge density will evolve like the left half side of Figs. [T^i 
and b: the charge density will be a narrow, positive function of x~ of width L, independent 
of x + , just as in the zero-temperature result (I3.47p . But then, at these times, 

QB(x) = (d t - V 2 )j°(x-) = -d^{l + d^)f{x-) (4.17) 

is the derivative <9_ of a function that is localized in x_, and so Q(x~) is a localized function 
whose integral vanishes. When resolving down to the scale L, the picture of fig. [SJo for j°(x) 
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> x 3 




> x 3 



FIG. 9: The space-time distribution of (a) Q@(x) = 9(Se(x) = (dt — V 2 )j°(x) and (b) its time 
integral Hq(x). Picture (a) is like fig. [3b but resolved down to the distance scale L. is negative 
in the gold region (marked — ) and positive in the other shaded regions (marked +). 



therefore becomes fig. [9^ for Q(x). The canceling positive and negative edges at early times 
blur together and disappear when we only resolve scales large compared to L. The time 
integral of fig. [9k, which defines £e as in (14.131) . is shown in fig. [9b- 

When we return to resolution scales large compared to L, where the positive and negative 
regions cancel at early times, then we approximate the x~ dependence of Q(x) by 5l{x~) as 
in our final result (II. lip . The coefficient of that S function will be the integral of Q(x) over 
x~ . Since S e was denned as the time integral of Q(x), this approximation is then just 



QQ(x) ~5 L (x-)£ (t=oo,x 3 ) 



(4.18) 



Our goal in what follows will be to use (I4.16P to compute Se(x) at t = oo, corresponding 
to the very top of fig. [9b . 



C. WKB approximations to Q\ 



Before we can evaluate the 4-momentum integral (14.16bj) that gives A(x, u), we first need 
formulas for the transverse bulk-to-boundary propagator Q R (q, u). This propagator is the 
solution to the linearized classical 5-dimensional equation of motion V/-F /J - = 0, which is 



d 2 u + ^rd u - 



fq 2 ~ ^ 



It will be useful to rewrite this equation as 



ft 2 + F) _ 



up 



Q±(uj, q,u) = 0. 



G±(u, q, u) = 0, 



(4.19) 



(4.20) 



where (as previously) q 2 = rj^q^qu = —u 2 + q 2 . In the high-energy limit, this classical 
equation can be solved using methods analogous to the semi-classical (WKB) approximation 
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in quantum mechanics, as discussed in Refs. |25|, |27| and in particular for the light-like case 
of q 2 = by Caron-Huot et al. [28]. 19 Here, we will need to examine small non-zero q 2 (with 
\q 2 \ <C uj 2 ~ E 2 ). In order to carefully understand the various scales at which different 
approximations are valid, we will go through the WKB approximation from the beginning. 
Treat u ~ k and substitute 

g ± = e <s-t+s +-) ( 421 ) 

in fl4.20p . where the exponent has been expanded formally in powers of 1/u for fixed u (with 
S n of order u~ n ). This gives 



7R 



C(q) 



u 



u 2 q 2 



1/4 



JS{q,u) 



(4.22) 



where we now use 5* as a short-hand notation for S-±, given by 

T 2ll/2 



du' 



[U' 2 q 2 



Q 



u: 



1/2 



/(«') 



(4.23) 



C(q) is an overall normalization factor not determined by the equation of motion. The 
approximation (" 14 . 2 2 f) is valid when the exponent S has large magnitude. 

We've written the answer in a form that's convenient for the time-like case q 2 < 0, which 
will be the most important later on. The choice of retarded propagator corresponds to 
taking the positive sign on the square root in (14.231) in this case. For space-like momenta 
q 2 > 0, the e %s analytically continues to an exponential suppression factor e - ' 5 '. 

Useful approximations to the integral (I4.23P for S depend on whether or not u is small 
enough that u' 2 q 2 can be treated as a perturbation in [u' 2 q 2 — q 2 ] 1 ^ 2 . The scale that 
separates different qualitative behaviors of S is therefore 



it* ~ 




(4.24) 



We discuss various expansions of f 14 . 2 3 1) in Appendix O The important limits for the present 
discussion are that 

2{-uq 2 ) 1 ' 2 



S 



for use when ti<t<l, and 



u 2 u 2 1 



\ + 0{u 2 ) + 0{^-) 



(4.25) 



S = ur(u) + lao- 1/2 ( 



(4.26) 



for use in the case u ^> u±. The constant c is defined by (I1.13p . lot(u) is defined as the 
result for S for the case o;=|q|, which gives 



Tanh" 1 v/w - Tan" 1 ^=-\n ( 1 + ^jj ) _ Tan" 1 (4.27) 

2 \l-yfuj 



19 Ref. [28| present their solutions in terms of the electric field E± instead of Q± oc A± . In our transverse- 
translation invariant problem, the relation is E± — iu>A± . There is also an additional difference in overall 
normalization: they do not normalize their solutions on the boundary like Q±. 
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1 I '- 1 

FIG. 10: The WKB approximation to G±{q,u) is valid in the shaded region parametrically far 
below the solid curve. Xhe solid curve corresponds to the matching scale tt ma tch 

discussed in the 

text, while the dashed line is the scale it* of (|4.24p . Ai and indicate the type of solution valid 
in the matching region for |g+| -C E^ 1 /' 3 and \q + \ S> E^ 1 ^ 3 respectively. 



with limiting cases 

( \ ft" 372 , U < I) (A rin\ 

r(u) ~ \ 3 1 4.28 
I — | m(l — u), u — y 1. 

The overall normalization C(q) of (I4.22p is fixed by the boundary condition Q±(q, 0) = 1, 
which lies outside of the region of validity for (I4.22p since S(q,0) = 0. Determining C(q) 
requires matching (I4.22p to a small- w solution of the equation of motion f)4.20p . We find 
C(q) in two complementary limits \q + \ <C E^ 1 ^ 3 and \q + \ ^> E^ 1 ^ 3 , which will be adequate 
for deriving our final result ( II. lip . We will discuss the appropriate scale for matching in 
each case below, but the result is summarized in fig. (TDJ The WKB approximation f !4.22p is 
a good approximation in the shaded region, parametrically far below the solid curve, which 
we will call w match (<?+)• 



1. £-V3 e 

The case E^ 1 / 3 <C \q + \ <C E (corresponding to E 2 ! 3 <C \q 2 \ ~ 4E|g + | <C E 2 ) will be the 
most important for studying the vast majority of charge deposition, which (as previewed in 
the introduction) will turn out to take place at distance scales x + <C E 1 ! 3 . 

For a given q + , let u ^> w m atch define the region of u where | ^S* | is large and so the WKB 
result (I4.22p is applicable. In order to determine C(q), we need to solve the equation of 
motion in the region u ~ Mmatch where the WKB approximation is marginal, which is the 
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region where \S\ ~ 1. If Mmatch <C «*, then eq. (I4.25P will give 

11 . 

Irl E\q+\ 

which will be consistent with (I4.24p for u* precisely when \q + \ ^> E~ x l 3 . 

Because m C C 1 where we need to do the matching, we can approximate 

q 2 - u 2 q 2 ~ q 2 (4.30) 

in the equation of motion (I4.20p . The complete approximation in this region is 



a 21 
d 2 -^- 

U 



g ± (cu,q,u)~0. (4.31) 



The boundary-normalized solution is just the vacuum solution (I3.16p . which we will write 
here in the form 

Of ~ iix^-uq 2 H^^-Auq 2 ) (u < <). (4.32) 

Using the asymptotic formula for the Hankel function to match to the WKB formula (I4.22p 
in the range w ma tch €«<it* where both are valid determines 



C(q) ~ e- i7r/ \-7iq 2 ) 1/2 . (4.33) 



2. \q+\ <C-ET 1/3 



When |g + | <C E 1 ^ 3 , we will see that \S\ ~ 1 at u matc h ^> u+, and so we turn to eq. (I4.26P 
for the WKB exponent S. Note that the second term in (I4.26P is of order 

w -l/2|g2|3/4 _ 

i? 1 / 4 ^! 3 / 4 <C 1 when |g + | <C E^ 1 / 3 , and so it can be ignored, leaving S — ujt{u). Then 
\S\ ~ 1 at 

Mmatch ~ ^ 2 / 3 ~ £~ 2 / 3 > it*. (4.34) 

This is consistent with (14.241) precisely when \q + \ <C E^ 1 ^ 3 . 

Because M ma tch 3> u± where we need to do the matching, we can approximate 

g 2 - u 2 q 2 ~ -u 2 co 2 (4.35) 

in the equation of motion (14.201) . Since also M ma tch <C 1, the complete approximation in this 
region is 

[d 2 u + uu 2 ] g ± {u,q,u) ~0. (4.36) 
The retarded, boundary-normalized solution is 

p> Ai(e" i7r / 3 Mw 2 / 3 ) , . . 

A.(0) (,I<<1 »- (4 - 37 > 

where Ai is the Airy function and Ai(0) = 3 _2 / 3 /T(|). 

One might worry that the analysis that led to (14.371) breaks down at u < u± <C w ma tch- 
But in that u range, g R is very close to its boundary value 1 and so the deviation from 
(I4.36P there will not affect the approximation (I4.37P at leading order in it* -C 1. 
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Using the asymptotic formula for the Airy function to match to the WKB formula (I4.22p 
in the range u ma , t ch < « < 1 where both are valid determines 



e «r/12 w l/3 



C ® * 2^AI(0) (438) 



and so 



e 



itt/12 



i , v ,> Ai{Q) u-^a;- 1 ^^ (« > Mmatch ). (4.39) 
Except for issues of overall normalization convention, this matching calculation is the same 



as the q = analysis of Ref. [28 
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D. Steepest descent analysis of A 



We now turn to using steepest descent methods to evaluate the integral ( I4.16b|) that gives 



A(x, u), just as we did at zero temperature in section [TH Dl As a qualitative preview of what 
we will find, fig. [11] is the finite-temperature version of fig. |HJ The horizontal axis is 

X + = x + -r(u). (4.40) 

Given our interest in distances x + ^> 1, eq. (I4.28[) for r(u) means that the difference between 
X + and x + is insignificant unless u is extremely close to the horizon. However, we will see 
that the behavior of A(x,u) as u— >1 is precisely what we want to get the large-time limit 
£e(£=oo,X3) that determines the charge deposition function via (I4.18p . We will see below 
that the x~ dependence of A is localized to x~ ~ — t(u) and so 

X + ~ x + + x' = 2x 3 . (4.41) 

There will be two different cases we will need to explore, corresponding to whether the 
saddle point of the q + integration probes the bulk-to-boundary propagator G R (q,u) above 
or below the curve u ~ u* of (14.241) and fig. [TIB 



1. Case A: Just like zero temperature 

At early times, we expect that the physics should be approximately the same as the 
zero-temperature case analyzed in section IIII Dl The zero-temperature bulk-to-boundary 
propagator (13.161) corresponds to the finite-temperature one when (i) u u+ and (ii) 
E^ 1 / 3 , so that (I4.32p applies. When these two conditions are satisfied, we may just take 
over the zero-temperature result (13.46bj) for A, 



A(x, e'°'- ^ Af (-|^; x -) O(^). (4.42) 



20 There is a typographic error in Eq. (A6) of Ref. [28| : The factor in big parenthesis in the left-hand 
equation should be raised to the 2/3 power. 
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x + - X ( u ) 




(EL 2 ) 1 ' 5 

FIG. 11: Qualitative picture of scales determining the behavior of A(x,u) at finite temperature. 
The horizontal axis is X + = x + — t(u), which is the same as x + except at the very bottom (u—>l) 
of the figure. The top curve u ~ x + /E indicates where the first wiggle is in A as a function of u, 
similar to the zero-temperature case of fig. The steepest descent approximation is valid in the 
shaded region (blue and magenta) below this curve (u 3> x + /E and X + -C E 1 / 3 ). In the lower 
shaded (magenta) region below the curve u ~ (x + ) 2 /EL, the field A is exponentially suppressed. 
The field is also exponentially suppressed in the region X + 3> E 1 ^ 3 to the right of the vertical 
black dashed line. 



From (13.401) . the saddle point is at 



uE 

|?+| = **~™. (4.43) 



Combining this with u± ~ ^/\q + \/E from (14.241) . the first condition u u+ is then 

which corresponds to the region above the solid curve in fig. [TTJ One may ignore the 
difference between X + = x + — t(u) and x + here because of r's relative insignificance away 
from the horizon. Using (I4.43p . the second condition \q + \ ^> E^ 1 ^ 3 will be satisfied if 

u > — and x + < E 1/3 . (4.45) 
E 

So the vacuum saddle-point result applies to the shaded region of fig. [U] that is above the 
solid curve and to the left of the black dashed line. 
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2. Case B: Falling into the black brane 

As we will see, the shaded region below the solid curve u ~ (a: + )~ 2 in fig. [11] will be 
determined by a saddle point with (i) u it*, and (ii) E~ l l z <C |g + | <C E. Combining 
(H23), p~26JL and P~33"j) with the fact that u 2 q 2 - q 2 ~ w 2 q 2 when it > it*, these two 
conditions give 

/ 2 \ 1/2 

r?R ~ „-^/4„-i/4 / \ J%,«) 



e- m/4 u- i/4 (-47rg+) 1/2 e ti+ w+We*«- n "\ (4.46a) 

where it) ~ g_ r(«) + SV(<?+; it) with 

it) ~ -q+ t(u) + §c£ 1/4 (-g + ) 3/4 . (4.46b) 

The e l<? - r ^) factor is important, and cannot be ignored, when u is very close to the horizon 
so that t(u) is large. When Fourier transforming from g_ to x~ , the effect of this factor 
will be to shift x~ by r(it). The high-energy approximation (" 13 . 18[) that we made in the 
zero-temperature case is modified to 

gf(q + ,q^u)^e^-^gf(q + ,u) (4.47) 

with 

g*t(q+,u) = e- i7T/4 u- 1/4 (-4nq + ) 1 / 2 e lS+ ( q+ > u) , (4.48a) 

giving 

A(x, it) ~ e iE[x ' +T{u)] J ^± Gf(q+, u) A[ 2) (q + ; x~+t{u)) e iq+x+ . (4.48b) 

Eq. f!4.48bp and the finite size L of the source region in x imply that A is localized to 

\x-+r(u)\<L, (4.49) 

and so x~ ~ — r(it) near the horizon. 

For a Gaussian source envelope (13. 20 p . approximating the integral (I4.48b|) by steepest 
descent requires extremizing 

S{q+, u) = -iS+(q+, u) - iq+x + + (q + L) 2 , (4.50) 

analogous to the zero-temperature case (13.391) . We will again treat L perturbatively and so 
find the extremum of 

S = -iS+ - iq + x + ~ -iq + X + - 2§c£ 1/4 (-g+) 3/4 , (4.51) 

which is at 

< 452 > 

One may now verify that the location of the saddle point satisfies the two requirements 
it 3> it* and E^ 1 ^ <C \q + \ C i? assumed for the propagator (I4.46P provided one is in the 
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Case B region of u ^> (X + ) 2 and 1 <C X + <C E 1 / 3 . Also note that the requirement 
l«So(g+)| > 1 is E/(X + f > 1, which is also satisfied when X + < E 1 / 3 . 

Expanding Sq to second order in small fluctuations q + —q+ about the saddle point, and 
then doing the Gaussian integral from (I4.48P and ( 14.46b|) . yields 

A(x, u) ~ e ^-H-M] e^V^M^) 1 / 2 A? } (£; aT+r(u)) f 2vr ^ 

To do a saddle point analysis, one should verify that there is a choice of integration contour 
that makes the neighborhood of the saddle point the dominant contribution to the integral. 
Having an explicit contour also helps one sort out exactly which branch one is on when 
evaluating the various roots in the derivation of (|4.53|) . We discuss the choice of integration 
contour in Appendix lEl 



3 - So (<?+,«) 



E. Final result for x 3 < E 1 / 3 



We are now ready to assemble our final result (11.111) for charge deposition up to distances 
of order E l l 3 . (The exponential tail at larger distances will be discussed in section HVFI ) As 
in the zero-temperature case, the derivative on A in eq. fl4.16aj) for Se(x) will be dominated 
by the term that hits e tEx , and so 



E e (x) 



A/" 2 
2E^ 



1 du 

«7 



\A(x, u) 



(4.54) 



analogous to (13. 22ft . We are interested in this result for t — > oo, which means arbitrarily 
large x~ and x + . Because of the localization (14.491) of A(x, u), non-negligible contributions at 
large x~ can only come from the near-horizon part of the u integration in (I4.54j) . where t(u) 
is large. Very large t(u) represents an exponentially-small region of u, and its contribution 
to the integral would be negligible if not for the 1/f factor in (I4.54[) . 
For t — > oo and «—)■!, our result (I4.53j) for A becomes 



8c 8 E 2 

Jx+y 



A 



(2) 



c 4 E 



(X 



;x +t(u, 



(4.55) 



Using (I4.28p to rewrite du/uf ~ dr in the u— >1 limit, and also using (14 . 4 1 [) in the same 
limit, 



£e(i=oo, x) 



Af 2 16c 8 E 3 
(2^a) 9 



POO 

/ dr 


A?( 


10 





c 4 E 



(2^ 3 ) 



-:x +t(u, 



(4.56) 



Shifting integration variable from r to x~ + r, and using the fact that \x~ + r\ is localized 
to L <C in the limit >oo of interest, we can use the definition (I1.17P and the 

result (I3.49P for the total charge Q per unit transverse area produced by the source (see also 
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Img + 



it/3/ *, 



Req + 



FIG. 12: A qualitative plot of the pole positions of Q± in the complex q + plane. The dashed line 
is proportional to e* 27r / 3 . 



Appendix IA 2j) to rewrite the last equation as 21 

£e(£=oo, x] 



- Uc^EL) 2 ( c 4 EL 
2Q v - J iff ( - 



(4.57) 



Using the relationship (I4.18P with the charge deposition function Q(x) finally produces 
our result fll.lip . As a check of the calculation, one may verify that this result satisfies 
J dt dx 3 Q(x) — \j dx~ dx + Q(x) = 1, as it should given the definition (11.9p of Q(x). 



F. The exponential tail of @(x) 

1. Relation to poles of 



We will now discuss the exponential fall-off of our result fll.lip for charge deposition for 
x+ > £ 1/3 . This requires evaluating A near the horizon for X + ^> E 1 ^ 3 , corresponding to 
the bottom of the white region in fig. [TTJ Near the horizon, we can always use the WKB 
formulas for C/J (see fig. ITU]) , and so we can use the integral expression (I4.48P for A, which 
we now find convenient to rewrite as 



A(x, u) ~ 
where Q± is defined by 



iE[x +t(u) 



gf(q + ,u)A^{q + ;x-+r(u))e^ 



2tt 



,u 



x+ 



t{u) 



U 



(4.5* 



(4.59) 



Attempting to evaluate this integral for large X + with saddle point methods will fail. But 
we can instead use the fact that the large-argument (X + ) behavior of a Fourier transform 
is determined by the singularities of that function in the complex plane (?+). 

In particular, the singularities of the bulk-to-boundary propagator are poles correspond- 
ing to quasi-normal modes of a vector field in the AdS-Schwarzschild background [25j, [29] , 
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Given our definitions (I1.18[) and (|3.19[) of A^ > (q+, <?-) and A { j^ > (q + ;x ), then J dx \A { j^ > (q + ;x 



(2) 



,(2)/ 
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shown qualitatively in fig. [12] for Q±. Remember that our convention is that q + = \((f — q°), 
and so the lower-half frequency plane (where the singularities of a retarded propagator should 
be) corresponds to the upper-half q + plane. 

We will review the origin of these poles in a moment, but first let's examine the conse- 
quence. For X + > 0, we can close the q + integration contour in (14.581) in the upper-half 
plane, giving a sum of residues from the poles q+=q+ \ each exponentially suppressed by a 

factor of e~ Im ^+ The dominant contribution at large X + will be from the pole q+=q^ 
closest to the real axis, giving 

-Im(q ( + ) )X+ 



A(x, u) 



le 



iE\x~ 



r(u)1 Res 



( 2 )^(!). 



(4.60) 



q+ will turn out to be small in the high-energy limit, so that it may be replaced by zero in 
the evaluation of the envelope function. Then 



\A\ 2 ~ Resffi( q { +\u)]A?(0;x-+T(u)) 
The same steps that we took from (14.551) to (14.571) then give 



-2Im(g(_ 1) )X+ 



E e (t=oo,cc) ~ AQL 2 



Res^(gf ) ,l)"e- 2Im ^)( 2 -3)^ (0) . 



(4.61) 



(4.62) 



For comparison, note that the zero temperature formula (13.161) for Q± has a branch 
point singularity at q+=0. This is why the zero-temperature result (I3.46bl) for A(x,u) falls 
algebraically rather than exponentially at large x + 22 



2. Scaling of poles with energy 

Poles in G±(q, u) occur when q is such that the normalization condition G±(q, 0) = 1 at 
u=0 causes G R (q, u) to be infinite for all other u. Turning this around, poles occur for q 
where finite solutions Aj_ (q, u) to the equation of motion with retarded boundary conditions 
at the horizon vanish at the boundary, Aj_(q, 0) = 0. For u ~ 1, we can use WKB methods 
to investigate solutions, but WKB breaks down and requires matching as u — > 0. When will 
this matching solution for u 1 cause A±(q, 0) = 0? 

For m < 1, the equation of motion (14.201) becomes 



■ 4Eq+-u 2 E' 
u 



A ± ~ (4.63) 



in the high energy limit. By the WKB analysis, the behavior of the retarded solution for 
u 3> M match is proportional to e lET( * u \ which for u 1 is 

S uV2E . (4.64) 



22 On the gauge theory side, the origin of singularities at q 2 = is the presence of long-lived massless 
excitations. At finite temperature, however, excitations generally have finite life-times due to interactions 
with the plasma. Only long-wavelength hydrodynamic excitations can have arbitrarily long life times. 
However, these will only couple indirectly to the large-momentum modes being considered here. On the 
gravity side, the coupling would be through loops, which are 1/-/V 2 suppressed, as in the discussion of 



long-time hydrodynamic tails of Refs. 22, 
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It will be simpler to analyze the question of when A(q, 0) vanishes if we can look at purely 
real solutions rather than complex ones. To this end, change variables to U = e~ m / 3 E 2 ^ 3 u 
so that the asymptotic behavior is 



A -2^3/2 



and the equation of motion is 



where 



A 



0. 



(4.65) 



(4.66) 



(4.67) 



This solution will have A±(q,0) = when the Schrodinger-like equation (14.661) has a zero- 
energy bound state solution that vanishes at the origin. That can happen for real positive 
a, which we label ai, 02, • • • starting from the smallest value that works. The corresponding 
pole locations are then 

Solving (14.661) numerically to find a\, we obtain a\ ~ 2. 141. 23 Correspondingly, the 
exponential decay factor in (I4.62p is 



e -2Im(g^ ) )(2 £C3 ) = g - ci(2x 3 ) / E 1 > ' 3 



(4.69) 



where 



a 1 



0.927. 



(4.70) 



Given the analytic structure of fig. [T21 readers may wonder what has become of the 
cut associated with the E l / A {— g+) 3//4 term in the WKB formula exponent (I4.46bl) . The 
formula for any contribution to a WKB exponent can only be trusted if its magnitude is 
large compared to 1, which for i? 1//4 (— g + ) 3//4 means when \q + \ ^> E^ 1 ^ 3 . But that means 
that \q + \ is large compared to the separation between poles in fig. [121 an d so the dense line 
of poles can approximate a cut. 



3. The residue 



Unlike the pole position, the residue will depend on u, and we are interested in the 
value near the horizon. From (I4.39p . we know that in this case there is a prefactor of order 
E~ l l & when \q + \ E^ 1 ^ 3 . Parametrically, the prefactor should be of the same order when 
\Q+\ ~ 1(2+1 ~ E~ l l 3 . So the behavior of the propagator near the pole should scale as 



£-1/6 



a — ai 



(4.71) 



23 We have double checked our analysis by also calculating the full, un- approximated bulk-to-boundary 
propagator numerically and verifying that we get the same scaling and pole locations at large ui. We did 



it by brute force, and one could likely find the poles more efficiently using the method of Ref. 29 1. 
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which we will use (I4.67P to write as 



1^1 - (4-72) 

\q+ - 9+1 

for some constant c^. By numeric evaluation of the full propagator ( 14. 59ft for smaller and 
smaller values of \q + — q+\ and u closer and closer to 1, we find 

c 2 ~3.2. (4.73) 

Using the residue from (I4.72p in the result (14.621) for E e (t=oo, x) then yields 

S e (t=oo, x) ~ 4Q^2 vj/( ) e ~^)/^ . (474) 
E 

Using (14.181) . this gives the x + 3> E 1 ^ 3 case of our final result ( II. lip . 
G. Revisiting the small-Q expression for G^^(Q-,u) 

We now return to discuss in hindsight the terms that we dropped in the small-Q form of 
Gaa(Qi u ) when replacing 

g*(u, k) rf v a v ~ ^-tj (too + ka 3 ) - — ^(1 - u) 1 ^ ' 2 {ka + coa 3 ), (4.75) 

Q*{u, k) rt v a v ~ -—^—(ia + ka 3 ) + — ^(1 - u)-^ 2 (ka + coa 3 ), (4.76) 

by 

g*(cu, k) vTa v -> -^(w + feo 3 ), (4.77) 

*) rTau (4.78) 

in section IIV A 11 We have seen in the transition from (I4.54p to (I4.57P how the 1/f factor 
was crucial to get a non-negligible result from near-horizon contributions, which in turn were 
crucial to get something for Eo(t=oo, x). That 1/f came from a factor g m of the inverse 
metric in our evaluation (14.61) of the 3-point function. It is only present for g 00 , not for g 33 . 
As a result, the Qf^ of ( 14.761) does not produce anything significant and can be dropped as 
in (14.781) . The terms of Qq" are multiplied by the 1/f factor in g 00 , but the second term in 
(I4.75P contains an explicit factor of 1—u, which cancels the near-horizon enhancement. So 
it too can be dropped, which is how we arrive at (14.771) . 



V. CONCLUSION 



We have shown how to use gauge-gravity duality for retarded 3-point correlators to solve 
a well-formulated gauge theory problem for studying the stopping of high-energy jets in 
strongly-coupled A^=4 super Yang Mills theory. Focusing on jets that carry R charge, we 
found more than one scale associated with the stopping of that charge, as described in the 
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introduction. It would be interesting to check whether our conclusions depend on the details 
of what we chose to study. In future work, one could study different observables, such as 
energy rather than R charge, and/or different sources, such as external gravitational fields 
rather than external R-charge fields. It would also be interesting to consider a source that 
is localized in the transverse direction and so look at the transverse spreading of the jet. 

Our final result (11. lip for charge deposition is exponentially suppressed at early times, 
before the first stopping scale. In giving (II. lip , we ignored details on distance scales < L. 
At early times there is a dipole contribution ~ 5' L (x~), as depicted in fig. [9j A moving, 
time-dependent dipole source might possibly produce a response that is not exponentially 
suppressed. In our approximations, the early time behavior was just the vacuum propa- 
gation of the excitation, which does not produce any charge diffusion. But there may be 
parametrically small corrections, dropped in approximations like (14.101) . that might produce 
dipole sources that do produce a small amount of diffusion originating from early times — 
that is, that produce effects which are suppressed but not exponentially suppressed. This is 
another possibility for further study. 
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Appendix A: Current conservation and total charge 

In this appendix, we will review how the Ward identity implies that currents are conserved 
outside of the source region. (As discussed in the main text, the current anomaly is not 
relevant to the particular field theory problem we have set up, and so we will ignore the 
anomaly in what follows.) We will also use the Ward identity to give an independent 
calculation of the total charge created by the source. Finally, we will show in detail how 
the Ward identity is respected by some of our main formulas and approximations in the 
gravity calculation. (We found such derivations very useful in the early stages of research 
debugging tool for our calculations.) 



1. General 



Start from the basic formula (12. 12[) for the response in terms of the 3-point function. To 
investigate <9 M (j 1 - 3 ^), we need to know Q^G^ ± . The Ward identity tells us that 

d tl G i ^ c)a ^{x 1 ,x 2 ;x) = f abc \s^{x-x 1 )G^(x 1 -x 2 ) - 5 {A) (x - x 2 ) G^ (x 2 - a*)] (Al) 
or equivalently [using Gr(— x) = G\(x)] 



(Qi,Q2;Q) = f abc Gi a (Q 2 )-G a /(Q 1 



yCt/3 i 



(A2) 



where both sides are of course multiplied by a momentum-conserving 5^ (Qi+Q 2 +Q) ■ The 
form (1A2I) of the Ward identity is familiar except perhaps for the details of retarded versus 
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advanced prescriptions. These can be quickly deduced by starting from the imaginary-time 
Ward identity and then analytically continuing in frequency according to (I2.18p . 24 
Because of the 5 functions in (1A1I) . 

Mj {3)tl (x)) Acl = \ / ^i^^Gf'^^i.^S^^^Oi^dN (A3) 
will vanish for x outside of the source region. 



2. The total charge created 

The total charge created is given by 

QW = J d 4 xd,(j^(x)) Aci . (A4) 
Using fl2T2D and the Ward identity flA"2|) . 

d * (j {c) ^))a , = \f abC I K a (Qi) - G a /(Qx)]AZi(Qi) A% cl {Q 2 ) e'^e-^*. (A5) 

JqaQ.i 



Then 



Q (c) = ~\r hc [ [Gf(Qi) - G^(-Qi)]A a * A (Qi) AlJQ,). (A6) 

JQi 

Now specialize to 

A a aiCl (q) = e a A a ±iCl (q) (A7) 

to get 

Q (c) = -\f ahc [ [Gi(Qx) - ^(-QOl^QxMi^QO 

^Qi 

= -if*[ pTiQ^ATAQi)^^), (A8) 

^Qi 

where we have used the relations G A (—q) = G A (q)* and ImG A (g) = |p NR (g). Here p NR 
is the spectral density with non-relativistic (NR) sign convention, related to the standard 
relativistic sign convention by 

p NR (g) = sign(gV el (g). (A9) 

To extract p±, we need the Green function G±. Since the source momenta Q\ are very 
large (of order k and so with components ^> T), we may use the vacuum result for the Green 
function, which is 

G ± = —J- lim z- x d z g ± = -— =- lim duQ± (A10) 



24 Alternatively, the derivation of (|Al|) directly in real time follows by applying to i 2 G^ bc " >a ^ l ' L (x\, x 2 \x) — 
0(t - t 2 )9(t 2 - t 1 )([[j^(x)J b ^(x 2 )],j aa (x 1 )]) + 9(t - t!)6(ti - t 2 )([\j^(x),j aa (x l )],j b ^(x 2 )}), using 
do9(t — ti) = 6(t — U) and the operator identity d^j^ = (ignoring the anomaly), and using the equal- 
time current algebra commutation relations [j a0 (t, x) , j bl1 (t, y)} = if abc j cll (t,y)S^(x — y). 
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with (I3.16P for Q±. This gives 



2^s 2 G 



[\n(uq 2 ) + 2 lE \. 



(All) 



The advanced prescription is q 2 — > q 2 + zesign(g°), giving 25 

n(-q 2 ) 



p NR (g) = 2 Im GUq) 



9sg 



6(-q 2 ) sign(q°) 



(A12) 



Using this in ( 1A8 



Q (c) = f abc I Ql 6(-Ql) AZ*{Qi) ^±AQi)- 



■2g : 



SG 



(A13) 



Now use the explicit form (11.21) for the source and / +3 = 2i: 



Q 



(3) 



vrA/1 



Ql sign(wi) 0(-Q?) |A(Qi - fc)| 2 - |A(Q! + fc)| : 



#SG 

/ g 2 signf.' l )^(-OT)|Af() J .-A-)|^ 



^SG ■/ Qi 

Since /c is large and A(Qi — k) localizes Q\ to be near k, we may set sign(wi) 
Q\ ~ 4Eq + , giving 



Q(3) 



8irEJ\f 2 



e(-q + )\q + \\A(q-k)\' 



9SG -Jq 

Then, for a transverse-translational invariant source 



, . e(-q + )\q + \\A(q)\ 2 . 

ySG Jq 



(A14) 
+1 and 

(A15) 



Q 



(3) 



SvrEATj [2dq + dq 



die 



9(-q + )\q + \\^ 2 \q + ,q^\ 



(A16) 



in agreement with (I3.49p . 



3. Ward identity for ([2734]^ 



Here we will check that the basic formula (I2.34bj) for the 3-point function satisfies the 
Ward identity. To investigate Q^G\ ± ^ we need Q^G^AQi x 5 )- The latter is the response 
A a (Q,x 5 ) to a boundary perturbation that is A^(Q,0) = Q^. But this boundary per- 
turbation can be gauged away (Ai — > Aj — diX) while remaining in A 5 =0 gauge by the 
x 5 -independent transformation X(Q, x 5 ) = i. The response to zero boundary perturbation is 
Aj = 0. Gauge transforming back gives the response A^ = Qn, independent of a; 5 . Therefore 



Q fl Q^(Q,x 5 ) = Q a (A 5 =0 gauge). 



(A17) 



2.1 



A good check of overall sign is that the spectral density p H should be positive for positive frequency. 
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Then, using Q = —Qi — Q 2 , f!2.34bj) gives 



iQ**G^ = -— ^ J d(x 5 ) v ^g ±± gP°(Q lp Q 1(7 -Q 2p Q 2(T )g^Q 1 ,x 5 )g^Q 2 ,x 5 ), (A18) 
where g- 11 - = eig IJ ej. 

Now consider the transverse equation of motion V/-F /J - = 0. For transverse-translational 
invariant sources, we will only need the Green function for Qi_l=Q2_l = 0. Using the fact 
that the metric coefficients depends only on u, we can write the equation of motion as 



0~3; (V=g g 11 g IJ {d J A L - d x Aj)) 
= -±=d,(^-gg ±± g^d 5 A L ) -g^g^q^A, 



and so 



Using this in ( 1A18f) gives 



(A19) 
(A20) 



iWGl^ = f rf(x 5 )[a 5 (v^^V 5 95^(gi,x 5 )) gi(Q 2 ,x 5 ) 

- Gt(Qi, x 5 ) d 5 (y/=g g^g™d 5 <#(Q 2| x 5 ] 



^ j ^ 5 )9 5 (v^^V 5 ^(Qi^ 5 )5 5 ^(g 2 ,a: 5 )) 



-gg^g^g^Qux^d.g^x 5 )) 



boundary 



(A21) 



where the last step implicitly assumes that the integral is sufficiently convergent that there is 
no contribution from the horizon (or u— >oo in the zero temperature case). At the boundary, 



g(q,x 5 )^l 



and 



gg^g 5 %g( q , X ^^-G(q). 



9sg R 



(A22) 



(A23) 



[For instance, in the AdS metric (13. 2p . the last is the usual expression (lAlOj) .] So we recover 
the Ward identity 

iQ>*Gf Xfl = G^{Q 2 ) - G^(Qi). (A24) 



4. Current conservation of ()3.8|) 



As a final example, consider the zero temperature expression (13. 8 aft in terms of A(x,u). 
Starting from ( I3.8a[) . 

d, (j^x)) oc r — [(jTdpdvAyA - A*(v^d u Ao)} ■ (A25) 
Jo u 
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The equation of motion ViF 11 - = for A is 



rrd„d v A = -udiA. 

Combining the last two equations, 

d, {j»(x)) oc / du [-{dlA*)A + A*dlA] 
Jo 

° ( ** a\ 

du da ( A* d a A \ 

- M->0 

A(x,0)d a A*(x,0) -h.c. 



(A26) 



(A27) 



Since A is proportional to the source on the boundary, the first factor A(x, 0) vanishes 
outside the source region, verifying that current conservation holds there. 



Appendix B: Evaluation of I(s) 

The first thing to notice is that the definition 

Hs) = J ^Gl^e-^j™ (Bl) 

of J(s) gives zero for s < 0. The argument is basically to close the integration contour in the 
lower half complex k plane, where has no singularities. Technically, one has to be a little 
careful because the convergence factor e~ eK does not converge for — 37r/4 < arg(/c) < — 7r/4. 
This problem can be avoided by first deforming the integration contour to run, for example, 
from e~* 77r//8 oo to the origin to e~ tn ^ 8 oo. At that stage, the e lKS factor produces a convergent 

2 

integrand, and one may drop the the now superfluous e~ £K convergence factor. Then one 
can close the integration contour at infinity, 

For s > 0, it's possible to evaluate the integral (IBlj) defining J(s) directly, by various 
contour deformation arguments and series expansions of Bessel functions. However, there is 
a simpler way using the equation of motion 

(%-^)Sx = (B2) 
satisfied by Q± at zero temperature. In terms of k = uq 2 , this equation is 

G± = Kd 2 K g ± . (B3) 

We can use this to rewrite ( 1B1I) as 

Hs) = f ^K[dlgl{K)]e-^e^. (B4) 

Integrating by parts twice, 

I( s ) = J ^ (s 2 k + 2%s)e-^e iK \ (B5) 
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which can be rewritten 

I(s) = i{s 2 d s + 2s)I(s). (B6) 
This differential equation is trivial to solve, giving 

J( s )oci e ^ s . (B7) 



s 



All that remains is to fix the overall proportionality constant by evaluating the original 
integral for some convenient value of s. This can be done for small s by the saddle point 
method of section HHP II (taking L = 0), or it can be done by evaluating the integral (IB 1 1) 
for large s by changing integration variables from k to A = us and then expanding the 
integrand in powers of 1/s. Either way, one obtains the result (I3.33p . 

Appendix C: Small Q form of Q R 

The low-Q behavior of the vector bulk-to-boundary propagator has been analyzed previ- 
ously [l2[ but not put into exactly the form that we need. What has generally been presented 
are the derivatives B u Qq and d u Qf- , whereas in this paper we want Qq- and Qf- directly, in 
A§=0 gau ge. It is easy to integrate, however, and determine the constants of integration. 
From Ref. 12], 26 keeping only the leading-order terms and for the sake of notational brevity 



writing A a for our • (cu, k) if" a 



d u A ~ —^(1 - u)-^ 2 (ka + coa 3 ), (CI) 
too — k z 



and 



d u A 3 = -— d u A ~ — — — (ka + ua 3 ). 

kj zoo — k z 1 — u z 



00 



2{%oo - k 2 ) 



1 + U 



(ka + ooa 3 ). (C2) 



Integration gives 



and 



Ao ~ C (w, k) - —^(1 - uf-^ / 2 (ka + ooa 3 ), (C3) 
too — k z 

A 3 ~ C 3 {u, k) + —^(1 - u)-^ 2 (ka + ooa 3 ), (C4) 
ioj — k 

where the integral of the (1 — u)~ %w l 2 j (1 +u) term from (1C2N has been dropped because that 
term's integral is sub- leading in powers of oo and k 2 . The integration constants Co and C 3 
are constrained by the fact that Aq and A3 must satisfy the Ao equation of motion, which is 

d 2 u A Q -—(kA, + ooA 3 ) = U 1 (C5) 
uf 

and so kCo + ooC 3 = 0. They are also constrained by the boundary normalization that 
A M —> dp. These constraints give (14. 7c p . 



2(, 



See specifically eqs. 73-78 of Ref. [12 1 
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R 



Appendix D: The WKB exponent S for Q 

Separately expand in powers of v! the factors of f~ l and [u' 2 q 2 — q 2 ] 1 ^ 2 in the second form 
of the integrand in (I4.23[) . Integrating term by term then yields 



S = 2«V2 ( _ g2) l/ 2 Fl (1. _ 1 , 1; |. *V j ^ ^ (Dla) 



*V 2 
where 

/3, /?'; 7 ; *, y) = £ £ (Dlb) 

m=0 n=0 

is the Appell hypergeometric function of two variables. The first term in this expansion 

gives flQSp . 

Rewriting the m sum in (IDlb|) as a hypergeometric function F = 2 -Fi gives the expansion 

00 r? 4- --n 4- ^ ^1 

5 = I(-V) 1/2 V 1 2."+4."+4. ^ J (D2) 

The standard hypergeometric transformation 

Th)T(B-a) 1 
F(a, /?; 7; 2) = F 7^ F ( (-z)- a F(a, a+l- T ; a+l-/3; -) + (a^) (D3) 

gives 

F(-I,n+i;n + f;^) 

= H±IW-i- _s_ , 1 g 2 xWV /2 r(n+f)r(-n-j) /-g 2 V + i 

71+ I 1 2 ' 71 3 ' n +3 'uV' V-?V 27rV2 ^2g 2y ) ' 

(D4) 

with which we can rewrite (ID 2 1) as 



M 2 « 



3 

4^(^)'S r( " + ^-"-t)(^ ! )"- < D5 > 



Expand the hypergeometric function as 



q > u 2n 



0000 / X\ / 2 \ m 

1„3/2|„|V^V^ ^"2^ / g \ „2n 



2 



V V 1 2>m u 2n (D6) 

-o m -( n - m + |) V^ 2 *? 2 ' 



m=0 
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Then rewrite the n > m part of the n sum as a sum over r = n — m, and use 

00 (-1) / n 2 \ m 
MEM 1 ? =^7 = ^ (D7) 



m! 

m=0 



and 



to obtain 



2r /-u „,/l/ 2 

r=0 4 



^ 3/2 E^I= / du'^—- 2 =T{u) (D8) 
r t 1 Jo 1 — It 



/ 2 \ 3/4 2 



m! \ ■( n — m + 7 / \ ircr" 

m=l \ n=0 4 / 

The first terms of this expansion give ( I4.26p . 

We have presented the expansion f1D9[) . useful when u 3> u*, as a series of tricks starting 
from the complementary expansion (IDlj) useful when u <C w* -C 1. It is also possible to 
derive ( ID9j) directly from the integral (I4.23P by appropriate expansions for u ^> u±. We will 
not reproduce here the full derivation of (ID9I) from this starting point, but the origin of the 
first terms (I4.26[) is easy to explain. Rewrite (I4.23[) as 

Jo u'^f(u') J, u' 1/2 /M 

where \x is an arbitrary scale with u t C <C 11. Then approximate as 



u' 1 ' 2 J, u>V 2 f{u>) 

'-'nonanalytic "T / wit 1/9/./ ,\ 

Jo «' x j(u') 

= ^nonanalytic + |<?| t(u), (Dll) 

where 

Sj =r^ [ "' v - g2l ',: 2 - [ " ,Vl ' /2 (D12, 

Jo U ' 

will not be analytic in q 2 (because a naive expansion of the integrand in q 2 leads to integrals 
with u— >0 divergences). A simple way to evaluate Snonanaiytic is to evaluate its derivative 
with respect to q 2 and then integrate back: 

U^nonanalytic _ \ j ^ / [ u H H \ 



d(q 2 ) " 2 J u' 1 ' 2 



00 , , KV-? 2 ]- 1/2 

2 lo u' 1 ' 2 



= -2- 3 / 2 c|g|- 1 / 2 (-g 2 )- 1 / 4 , (D13) 

and so 

^analytic * |c| q | (- 1^)3/4_ (D14) 

Approximating |q| ~ u; in (IDllj) and (ID14[) then gives the leading terms shown in (I4.26|) . 
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Appendix E: Saddle point integration contours 



When making a saddle point approximation, one should choose an integration contour 
such that the contributions to the integral are negligible everywhere except in the neighbor- 
hood of the saddle point. Fig. [13] shows examples of integration contours that do the job 
for Case A and Case B analyzed in section HVDl (We consider here a Gaussian envelope 
function for the sake of concreteness.) The dashed line depicts the line of poles discussed in 
section [TV F| which also serve as the location for the cut of the (—uEq + Y^ 2 and E x l^[— g + ) 3//4 
terms in the WKB exponents (14.251) and (I4.26p . The integrand of the q + integral in (I4.48b[) 
will be exponentially suppressed compared to the saddle point (shown by the large dot) in 
the interior of the shaded region. The various circles indicate different scales for \q+\, as 
labeled. 

The reasons for the suppression are different in different regions. We will go through 
Case B as an example. We get exponential suppression when 

S(q+, U) ~ -if C £V4 ( _ g+) 3/4 _ %q+X+ + {q+L) 2 (E1) 

has a positive real part ^> 1. In the interior of the smallest circle (|g + | -C \q\\ ~ E/(X + ) A ) in 
fig. [13b, the — iE l l 4 (— q+) 3 ^ 4 term dominates, and the shaded part shows where its real part 
is positive. In the next annulus (E/(X + ) 4 <C \q + \ <C u 2 E), the — iq + X + term dominates, 
and the shaded part shows when its real part is positive. In the next annulus out (u 2 E <C 
\q + \ <C X + /L 2 ), \q + \ is large enough that the expansion (lEip is no longer appropriate, and 
we should switch from (14.261) to (I4.25p . giving 

S(q + ,u) ~ -iA{-uEq + yi 2 - iq + X + + (q + L) 2 (E2) 

as in (I3.39p . But the — iq + X + term still dominates. Finally, beyond the outermost circle 
(|g+| 3> X + /L 2 ), the (q+L) 2 term dominates. 
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